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HAUPTAUFSÄTZE 


Die Prandtl-Glauertsche Näherung als Grundlage für 
ein lterationsverfahren zur Berechnung kompressibler 


Unterschallströmungen. 
Von W. Hantzsche in Braunschweig. 
(Aus dem Institut für Gasdynamik der Luftfahrtforschungsanstalt Hermann Göring E. V., Braunschweig.) 


Für ein Iterationsverfahren zur Berechnung der kompressiblen Strömung, 
das von der Prandtlschen Näherung für dünne Profile als erstem Schritt 
ausgeht, wird die dritte Näherung berechnet. Unter Verwendung der kon- 
formen Abbildung, die das Außengebiet des Profils in der nach Prandtl 
affin verzerrten Strömungsebene in das eines Kreises überführt, läßt sich 
die dabei zu integrierende inhomogene Differentialgleichung allgemein lösen. 
Für die Strömung um die nichtangestellte Ellipse und den nichtangestellten 
Kreisbogen wird die dritte Näherung vollständig angegeben und numerisch 
ausgewerte. Der Kompressibilitätseinfluß steigt gegenüber der zweiten 
Näherung noch weiter om 


Die Prandtl-Glauertsche Näherung für die kompressible Potentialströmung besteht 
in der Annahme kleiner Zusatzgeschwindigkeiten zu einer Parallelströmung als Grundströmung. 
Die dadurch erzielte Linearisierung der Differentialgleichung des Potentials bzw. der Strom- 
funktion einer adiabatischen Gasströmung gibt den Kompressibilitätseinfluß für dünne und 
flach angestellte Profile gut wieder. Auch für dickere Profile stellt die Prandtl-Glauert- 
sche Linearisierung eine recht brauchbare erste Näherung dar. In einer früheren, gemeinsam 
mit H. Wendt verfaßten Arbeit‘) wurde nun für die ebene Strömung eine Methode der 
sukzessiven Approximation für die Stromfunktion angegeben, die die Prandtl-Glauertsche 
Näherung als ersten Schritt hat. Wirklich durchgeführt ist dort aber nur der zweite Schritt. 
Die vorliegende Arbeit gilt der Berechnung des dritten. 

Die Schwierigkeit bei der praktischen Durchführung von schrittweisen Näherungen be- 
steht neben der Erfüllung der Randbedingungen in der Angabe des allgemeinen Integrals der 
bei jedem Schritt neu zu lösenden inhomogenen linearen Differentialgleichung. Bei der Be- 
handlung des zweiten Schrittes konnten wir zeigen, daß sich die allgemeine Lösung der in- 
homogenen Gleichung explizit angeben läßt, so daß nur noch die praktisch natürlich oft auch 
schwierige Befriedigung der Randbedingungen durch eine geeignete homogene Lösung nötig 
blieb. Die Behandlung der Differentialgleichung des dritten Schrittes zeigt nun, daß man 


1) W. Hantzsche und H. Wendt: Der Kompressibilitätseinfluß für dünne wenig gekrümmte Profile bei Unter- 
schallgeschwindigkeit. Z. augew. Math. Mech. Bd. 22 (1942), 8. 72 (vgl. dort auch den Hinweis auf A. Busemann). 
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auch hier die allgemeine Lösung angeben kann, wobei sich allerdings im Gegensatz zur zweiten 
Näherung von den auftretenden Quadraturen eine allgemein nicht explizit durchführen läßt. 
Sie stellt aber praktisch in vielen Fällen ein Integral über eine rationale Funktion dar. 

Die dritte Näherung wird für die Schar nichtangestellter affın verzerrter Kreisbogen und 
die nichtangestellte Ellipse berechnet. 


8 1. Aufstellung und Transformation der Differentialgleichungen für die ersten 
drei Schritte. 

Zur Herleitung der für die einzelnen Näherungsschritte gültigen Differentialgleichungen 
hatten wir für die auf ein kartesisches x y-Koordinatensystem bezogene Stromfunktion y einer 
ebenen adiabatischen Gasströmung um ein Profil eine Reihe in einem Parameter Ä angesetzt, 
der die Abweichung der Profilform von der nichtangestellten Strecke charakterisiert?): 


Dabei liegt die Anströmung U in Richtung der negativen x-Achse. Der ‚Zusammenhang 
zwischen der Stromfunktion y und den Geschwindigkeitskomponenten u, v in ' Richtung der 
x- und y-Achse und der Dichte o wird durch die Beziehungen: 


0% 


vermittelt (0, Dichte im Unendlichen). 
Geht man mit dem Ansatz (1) in die Differentialgleichung 


u? 
(1- =) -5)=0 
(a Schallgeschwindigkeit) der Stromfunktion einer adiabatischen rotationsfreien Gasströmung 


ein, so ergeben sich durch Nullsetzen der Koeffizienten von A, 4? und 4? für die ersten drei 
Funktionen g (x, 4), h(x,y) und f(x,y) der Reihe (1) die Differentialgleichungen ®): 


guy — [x a +1] u’ 9x gu 
Mit a ist die Machsche Zahl der Anströmung bezeichnet | == Zn, a, Schallgeschwindigkeit 


im Unendlichen)), u.’ bedeutet das Verhältnis der spezifischen Wärmen bei konstantem 
Druck und Volumen. Zur Abkürzung ist weiter gesetzt: 
1 a 
3 a u? 


Führt man an Stelle von g(z, 4), h(x,y) und f(x, y) die Abschnitte der Stromfunktion y ein: 


T 

1 


2) Für welche Strömungen sich % wirklich nach einem Parameter A entwickeln läßt, kann allgemein nicht ent- 
schieden werden. Auf eine Verallgemeinerung von (1) werden wir in Beispiel B von $ 3 geführt. 
3) Vgl. S. 74 von 
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so überzeugt man sich leicht, daß das System der Differentialgleichungen (2) mit dem folgen- 
den wird: 


\ 


| (3e). 


U xy 


Das Näherungsverfahren wird nun in der Weise durchgeführt, daß wir drei im allgemeinen 
von y®, 9 verschiedene „Stromfunktionen* (z, 9,4), (z,9,4), 
bestimmen, die den Randbedingungen von y genügen: 


a) (x, y,A)—=const auf dem Profil, 


b) IT md, 


c) _ Abfluß an der Hinterkante des Profils bzw. Staupunkt 
und die Differentialgleichungen (3) erfüllen: 


—= — U im Unendlichen (j=1,2,3), 


(2) 


(1)\2 
3 Yx 24 1 


Es ist zu erwarten, daß die Näherungen y“”, w“ bzw. y“® in ihrer Entwicklung nach A mit 
der Reihe (1) für y bis zum linearen, quadratischen bzw. kubischen Glied übereinstimmen, 
da wir in den Differentialgleichungen (4b) bzw. (4c) gegenüber (3b) bzw. (3c) auf der rechten 
Seite nur Fehler in jeweils höheren Gliedern von A machen '). 

Die Näherungsfunktion y“ (x, y, A) ist offenbar die Prandtl-Busemannsche Näherung 
für die Stromfunktion. Die. Differentialgleichung (4b) wurde früher behandelt, wir wollen 
jetzt auch noch (4c) integrieren. 

Dazu transformieren wir die Differentialgleichungen (4a, b, c) auf geeignete unabhängige 
Variable. Durch die affine Verzerrung 


erreicht man, daß auf der linken Seite der Differentialgleichungen der Laplacesche Eat 
auftritt. Das Profil der €n-Ebene geht aus dem der &y-Ebene durch diese Verzerrung her- 
vor. Das Außengebiet des verzerrten Profils werde nun konform auf das des Kreises 


4) Für den Fall überall gleicher Zusatzgeschwindigkeiten, also staupunktsloser Strömung‘ kann dies bei einigen 
in praktischen Beispielen im allgemeinen erfüllten Voraussetzungen über die Profilformen bis zur dritten Näherung 
streug bewiesen werden. 
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r=yr’+y’=a'einer r,9-Ebene abgebildet. Bezeichnen wir mit £ bzw. 3 die komplexe 
Variable &+in bzw. g-+iYy, so sei die Abbildungsfunktion durch die Beziehung 


gegeben. In der Umgebung von 3=x sei diese Abbildung die identische. Mit 
nehmen dann die Differentialgleichungen (4a), (4b), (4c) die Gestalt an: 
0 
AK u? . . . . . . . (6b), 
U 
xD x” 
(6e). 


Für die Funktionen (r,9,4), K®(r,9,4), K“(x,9,4) haben wir offenbar die Rand- 
bedingungen: 
a) KW)=const auf dem Kreise r=a, 


b) 0, U im Unendlichen (j=1,2,3), 
e) Ti 0, Tu 0 im Bildpunkte der Hinterkante des Profils. 


$ 2. Allgemeine Integration der Differentialgleichung des dritten Schrittes. 
Für die Integration der Gl. (6c) brauchen wir zur Berechnung der rechten Seite die 
Struktur der Funktionen X‘ (£,n) und X” (£,n) (die Abhängigkeit dieser Funktionen vom Para- 
meter A werde hier und im folgenden nicht mehr ausdrücklich angegeben) und ihrer ersten 
partiellen Ableitungen nach &£ und „7. Wir geben diese daher in den Abschnitten A und B 
an. Unter C wird die rechte Seite von (6c) aufgestellt. Der Absatz D bringt die Integration. 


A. Die Gestalt der Funktionen X’ (£,n) und (£,»). Bezeichnen wir mit das kom- 
plexe Potential der inkompressiblen Strömung in der ry-Ebene um den Kreis r=a bei der 
Anströmgeschwindigkeit — U und den Bildpunkt der Hinterkante als Staupunkt, so ist K“” 


offenbar: 


Mit $m (bzw. Re) wird der Imaginärteil (bzw. Realteil) einer komplexen Größe bezeichnet. 
Wie in der in Fußnote ') zitierten Arbeit gezeigt wurde, und aus der Berechnung von K“” 
nochmals mit hervorgehen wird, hat K“ dann die Gestalt: 


Dabei bedeuten 


. (8a) 
und weiter ist zur Abkürzung gesetzt: 
2 
. (8b). 


de: 
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Der Strich gibt den Übergang zur konjugiert komplexen Größe an und die Funktion Qß@) 
bestimmt sich durch die Randbedingungen von K“®. 
Die Ausdrücke für X“ (&,n) und X (£,n) ergeben sich dann wegen (5) aus (7) und (8): 


[637 


X 
U = 


B. Die Berechnung der ersten und zweiten partiellen Ableitungen von X‘ (£, 7) und X (£, n). 
Zur Berechnung der Differentialquotienten von X‘ (&,n) und X(£,n) nach & und n wollen 
wir zunächst einige Hilfsbetrachtungen über das Rechnen mit komplexen Ausdrücken hinter 
dem Zeichen Re und Im anstellen. 


Für eine analytische Funktion hat man die Beziehungen: 


Durch Übergang zu konjugiert imaginären Größen bekommt man daraus: 


‚df@- 
dh 


Geht man in den analytischen Funktionen f(;) un 
FO, so. ergibt sich leicht: 


— — — 


(3) mit {={£(4) zu der Variablen £ über: 


F dF 
dF 
- 
Damit erhält man die Relationen 
lfd dfd 
(10) 
dfd . dfd 
und 
9 _m,[dfdi, ‚dAdz 
Rechenformeln, deren allgemeinen Aufbau man leicht übersieht. 
Die Ableitungen von X“’(&,n) ergeben sich aus (9) nun sofort mit (10) zu: 
P=Imm, | 
ERROR 
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setzen. Die Ableitungen von X‘“”(£,n) erhalten wir aus (9) und (11): 


17 = Im N], 
IH MM 1, 
| 
am 
_ IM | 
na 
. (136). 
dad} 


C. Die Aufstellung der. rechten Seite von (6c). Da für die komplexen Größen ® und Q 
27m 


ist, können die einzelnen Glieder der rechten Seite von (6c) in der folgenden Weise ge- 
schrieben werden: 


_ MR OMT 
ax) x? 


[MRHMN], 


EM 


x® 
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x» xD 
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x» 2 


«3» 


Damit bekomnit man aber für 4K 


Ü die Gestalt: 
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+, ENN +3, 
+2, 


wo die v; j=1,...,11) folgende Abkürzungen sind: 


4 2 


. (16). 
1 1 2 1 2 


1 


1 


D. Die Integration von Gl. (6c). Zur Integration von (15) führen wir als unabhängige 
Variable 3 und 3 ein’). Wegen 
K“ 0° K (3) 0° 
kann das allgemeine Integral von (15) angegeben werden, wenn die Doppelintegrale in 3 
und 5 über die einzelnen Terme der rechten Seite bekannt sind. Diese lassen sich aber 


weitgehend ausführen oder sind zumindestens auf einfache Quadraturen zurückführbar, weil 


-5- a ist und die einzelnen Glieder der rechten Seite von (15) ebenfalls in ein Pro- 


dukt einer Funktion von 3 und einer in 3 zerfallen. So hat man z. B. 


er 
dt 
Entsprechend ergibt sich: 
_dt TE ‚(d& AM 


5) Die Rechtfertigung des formalen Rechnens mit 3 und 5 als unabhängige Variable findet sich auf S. 77 von !) 
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dad da 


Schreibt man nun K“ selbst an und ordnet die Terme, so bekommt man leicht: 


+, MMS+r, ZW +, FSEM+r (17). 
+r,8P +, WIN + MM 


Bis auf ein Integral haben sich also die noch nieht ausgeführten Quadraturen herausgehoben. 
F(3) bestimmt sich dann aus den Randbedingungen für K‘”. Die untere Grenze des auf- 
tretenden Integrals kann beliebig gewählt werden. Eine Änderung der Integrationskonstanten 
bedeutet das Hinzutreten einer Funktion in 3 und damit nur eine Änderung von F(3). Das 
Integral selbst wird sich praktisch oft auswerten lassen. Mit Rücksicht auf die allgemeine 
Gestalt von ® (3): 


a 
ist bei in 3 rationalem £(3) der Integrand nämlich stets eine rationale Funktion von 3. 


8 3. Beispiele. 

A. Die Berechnung der dritten Näherung für die Schar der affin verzerrten nichtangestellten 
Kreisbogen. Die Durchrechnung der einzelnen Näherungen setzt, wie der vorangehende 
Paragraph zeigt, die Kenntnis der inkompressiblen Strömung um das Profil der &7-Ebene und 
die Abbildung des Außengebietes dieses Profiles auf das eines Kreises der ry-Ebene voraus. 
Gibt man ein Profil in der Strömungsebene (x y-Ebene) vor, so wird es im allgemeinen 
schwer sein, die inkompressible Strömung um das affin verzerrte Profil und die zugehörige 
Abbildungsfunktion zu finden. Dazu kommt noch, daß die affine Verzerrung von der Mach- 
schen Zalıl abhängig ist, für jede Machsche Zahl die Rechenarbeit also neu durchgeführt 
werden muß. Es ist daher praktisch einfacher, sich das Profil der affin verzerrten &n-Ebene 
vorzugeben. Dann hat man zwar den Nachteil, daß die daraus gewonnenen Profile der 
x y-Ebene mit der Machschen Zahl sich ändern, kann aber in der &n-Ebene das Profil so 
wählen, daß die inkompressible Strömung und die Abbildungsfunktion gegeben sind. 

Als Beispiel für die Berechnung der dritten Näherung wollen wir zunächst in der 


£n-Ebene den nichtangesteliten Kreisbogen der Wölbung — wählen. In der x y-Ebene be- 


kommen wir dann eine Schar affin verzerrter Kreisbogen der Wölbung 7 und zwar ist für 
. a=1l(a=0) das Profil der «y-Ebene offenbar ein genauer Kreisbogen, mit wachsendem a 
erhält man Ellipsenbogen, die im Grenzfall a>1 in einen Parabelbogen übergehen. 

Die Abbildungsfunktion 
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bildet den Kreis r=a der 

3-Ebene _auf einen Kreis- 

a 7„-$ bogen AB der Pfeilhöhe f 
in der £n-Ebene ab (Bild 1). 

, A A und B sind die Bilder der 

- y zy-£bene Punkte A’ und B’ des Kreises 
der j-Ebene. In Y- 

‘bene erhalten wir dann 

einen Ellipsenbogen A,B, 

Pas) 4 mit der Pfeilhöhe fu. Die 


ild 1a. Bild 1e. = 
Wölbung ı dieses Bogens 


ist, wie man sich leicht überzeugt, mit dem Parameter f durch die Beziehung 


verknüpft. 


Als Entwicklungsparameter verwenden wir zunächst f, das Ergebnis läßt sich ja dann 
leicht auf die Wölbung 7 umschreiben. 

Die erste Näherung K“’ ist nach (7) durch das komplexe Potential UW(3) der inkom- 
pressiblen Strömung um den Kreis r=a der x,9-Ebene gegeben. Die Anströmgeschwindig- 
keit ist — U, die Staupunkte liegen bei A’ und B’: 


a? 
Die Gestalt von K“ ist nach (8) dann durch die Kenntnis von © (5), M(5) und Q (3) festge- 


legt. Dabei berechnen sich &(5) und M(3) nach (8a) aus der Abbildungsfunktion (18) und 
dem komplexen Potential (19): 


a 
(20) 


Die analytische Funktion Q (3) muß so konstruiert werden, daß K“ die Randbedingungen 
(a), (b), (e) von S. 188 erfüllt. Dies wurde in unserer ersten Arbeit unter der für die zweite 
Näherung zulässigen Beschränkung auf die quadratischen Glieder in f getan. Mit Rücksicht 
auf die Verwendung des Ausdruckes von K“’ in der dritten Näherung brauchen wir aber 
jetzt Q bis zu den Gliedern in f? genau: 


(21). 


Die zweite Näherung hat dann (bis zu den dritten Potenzen in f genau) die Gestalt: 
K? 1 a? In a? 
In — In —— 


= 1 +3 a? + a? | 3 
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Die Ermittlung von Q\ gelingt dabei leicht mit dem schon in der vorangegangenen Arbeit 
verwandten Konstruktionsprinzip. Eine in 5 und 3 analytische Funktion H (3,3) hat auf dem 


2 
Kreise r=a (33 offenbar dieselben Randwerte wie H =) oder H Der Imaginär- 


2 2 
teil H (3,3) kann deshalb entweder durch Jm oder H ) für r=a kom- 


pensiert werden. Damit gelingt es, die Randbedingung (a) zu erfüllen. Man muß dabei 
allerdings darauf achten, daß im Außengebiete des Kreises keine Singularitäten auf- 


treten bzw. diese (wie in unserem Beispiel Pole im Unendlichen durch Terme der Form 


> wieder beseitigen. Nach der Befriedigung der Randbedingung (a) 
bieten dann rd beiden anderen ker zT keine Schwierigkeiten. Der erste Bestand- 
; der zweite mit x, und x, behaftete hat im Unend- 
lichen ee Me Ableitungen. Der Staupunkt (Bedingung (c)) wird schließlich durch 
Addition eines Vielfachen von Jmi Int erfüllt. 


Zur Ermittlung von K“ brauchen wir [vgl. (17)] noch die Größen ®, N und S, die wir 
nach (13b) leicht berechnen können. Wir geben sie in ihrer Entwicklung nach f jeweils nur 
bis zu der Ordnung an, wie wir sie im Ausdruck von K“ benötigen, um alle Terme in f, 
f” und f* zu erhalten: 


3a’n, 
‚In a4... 


Das in (17) auftretende Integral u sich leicht ausführen. Bestimmt man weiter F(z3) nach 
dem eben bei der Berechnung von K“” erwähnten Konstruktionsprinzip,“ so erhält man 
schließlich für K“ den ren Ausdruck: 


33 

ln 


33 

In dieser Darstellung für K” tritt Rs für 3= +a ein Pol auf. Man überzeugt sich aber 
leicht, daß K“ an diesen Stellen endlich bleibt und dort auch stetige erste Ableitungen hat. 

Mit der Kenntnis von K“” ist nach (5) und (18) das Strömungsfeld in der x y-Ebene 
‚bekannt. 

Zunächst werde die Geschwindigkeit im höchsten Punkte des affın verzerrten Kreis- 
bogens (Maximalgeschwindigkeit) der Wölbung r angegeben. Der Zusammenhang von r und f 
wird durch (18a) geliefert. Man hat für die Maximalgeschwindigkeit bis zu Gliedern in +’ 
genau 
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58 44 10 1350 
Der Auftriebsbeiwert c, = > ergibt sich zu: 
: Tiefe 

13 25 1 1 


(35). 


25 7 
In Bild 2 ist die Maximalgeschwindigkeit für den Fall = als Funktion der Mach- 


schen Zahl aufgetragen. Das Verhältnis der spezifischen Wärmen wurde dabei 1,4 gewählt. 
Die im Geschwindigkeitsbild eingetragene Kurve € gibt an, wann die kritische Geschwindig- 
keit erreicht wird. Als erste, zweite bzw. dritte Näherung ist dabei jeweils die Entwicklung 
bis zum Gliede in z, 7? bzw. 1°? bezeichnet. 


Wie wir schon Seite 192 angaben, deformiert sich bei konstanter Wölbung 1-5 bei 


der Änderung der Machschen Zahl von 0 bis 1 der Kreisbogen über Ellipsenbögen zum 
Parabelbogen. Bis zum Erreichen der kritischen Machschen Zahl (a=0,575) ist diese 


Profildeformation aber sehr gering, wie man aus Zahlentafel 1 ersieht, wo * in Abhängigkeit 
von 7 aufgetragen ist. (x-Achse in der Sehne, y-Achse Symmetrieachse des Kreisbogens, t Tiefe.) 


Zahlentafeli. 


-U 
a=0 a=0,3 a= 0,575 
# t t t t 
/ 0 0,1 0,1 0,1 
0,1 0,09615 0,09614 0,09610 
0,2 0,08452 0,08448 0,08435 
4 / 0,3 0,06491 0,06483 0,06461 
0,4 0,03693 0,03685 0,03662 
\ 0,5 0 0 
- 
m Bild 2 (links). Maximalgeschwindigkeit am affin verzerrten Kreisbogen 


20 (Wölbung „) in Abhängigkeit der Machschen Zahl. Kurve 1,2,3 erste, 


02 zweite, dritte Näherung 1”, ..., ©’ Prandtl Glieder bis r,r?,r3, .. 
vollst, Lösung. 


Berechnet man sich einmal die Geschwindigkeitsverteilung bei inkompressibler Strömung 
für den bei der Machschen Zahl a sich einstellenden Ellipsenbogen bis zu den Gliedern 
in a’ genau, so erhält man: 


Man sieht daraus, daß im inkompressiblen Fall die Maximalgeschwindigkeit (£ 4 bis zur 


dritten Näherung von der Profildeformation nicht beeinflußt wird. Die Kurven 1, 2 und 3 
der Darstellung 2 dürften deshalb die Abhängigkeit der Maximalgeschwindigkeit von der 


Machschen Zahl für einen Kreisbogen der Wölbung io gut wiedergeben. 


Zum Vergleich mit der Prandtlschen Näherung ist in Bild 2 noch die nach der 
Prandtlschen Regel berechnete Maximalgeschwindigkeit am Kreisbogen nach r entwickelt 
aufgetragen. Es zeigt sich dabei, daß die dritte Näherung (Glieder bis einschließlich ı’) die 
vollständige Lösung schon gut annähert. 
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Bild 3. Auftriebsbeiwert am affin verzerrten Kreisbogen 

(wörbung 4) in Abhängigkeit der Machschen Zahl. 
Auftrieb 

U?2.Tiefe 


U Geschwindigkeit 


Schallgeschwindigkeit Unegdlichen. 


Bild 4. Druckverteilung am Kreisbogen der Wölbung 5; 
in erster bis dritter Näberung. 


Bild 5. Inkompressible Druckverteilung am Kreisbogen 
der Wölbung u in erster bis dritter Näherung. 


Bild 6. Druckverteilung am Kreisbogen der Wölbung = 
in dritter Näheruug. 


724 Bild 7. Druckverteilung für die inkompressible Strömung - 
-76- an den affin verzerrten’ Kreisbögen (M=0, M = 0,575) der 


1 
Bild 5. Wölbung ;;- 


In Bild 3 ist der Auftriebsbeiwert Fr nach (25) auch in Abhängigkeit von der Mach- 


schen Zahl a bei konstanter Wölbung n angegeben. Die in inkompressibler Strömung durch 


die Profildeformation eintretende Auftriebsänderung ist ebenfalls sehr gering: 
und wird in Bild 3 durch die gestrichelte Kurve dargestellt. 
In Bild 4 ist für die Machsche Zahl 0,575 (kritische Machsche Zahl der Anströmung) 


die Druckverteilung über der Sehne des nichtangestellten Kreisbogens der Wölbung ri ge- 


zeichnet. Aus dem entsprechenden Bild für die inkompressible Strömung um einen Kreis- 
bogen (Bild 5) sieht man, daß dort mit der dritten Näherung bereits die exakte Lösung in 
Zeichengenauigkeit angegeben wird. 

Bild 6 gibt schließlich noch in dritter Näherung die Druckverteilung für a=0, a=03 
und a=0,575 in einer Darstellung. Dazu ist allerdings wieder zu bemerken, daß für die 
drei Machschen Zahlen die Profile nicht ganz übereinstimmen. Ihre Abweichung ist aber, 
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2757 73 24 Z 48 (74 ARE -16 
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-16 
Bild 4. 
u 
- -/Näherung 


Ba. Nr Hantzsche, Die Prandtl-Glauertsche Näherung 197 


wie bereits erwähnt, gering. Bild 7 zeigt, daß die Profildeformation bei der inkompressiblen 
Druckverteilung (in dritter Näherung) nicht ins Gewicht fällt. Es wird dort der exakte 


Kreisbogen (@a=0) der Wölbung 2 mit dem der Machschen Zahl a =0,575 entsprechenden 
affın verzerrten gleicher Wölbung verglichen. 


B. Die Berechnung der dritten Näherung für die nichtangestellte Ellipse. Als weiteres 
Beispiel wollen wir jetzt die nichtangestellte Ellipse behandeln. Die am Anfang von A. er- 
wähnten Schwierigkeiten, daß sich das Profil im allgemeinen mit der Machschen Zahl 
deformiert, tritt hier nicht auf, weil eine affine Verzerrung eine Ellipse wieder in eine solche 
überführt. Die Funktion 


bildet den Kreis r=a der 3-Ebene in eine RE mit dem Achsenverhältnis der 


5 
£-Ebene ab. In der «y-Ebene bekommt man dann eine Ellipse mit dem Dickenverhältnis 
e=u a’ + . 


Die erste Näherung ist offenbar: 
a?’ 

Die Funktionen (3) und M(3) sind dann mit A=a?’ — b*: 


A 


A 
und Q (3) bekommt, wie wir schon in unserer früheren Arbeit zeigten, die Gestalt: 


mA) mA 


a? 
Zur Berechnung der dritten Näherung benötigen wir nun die leicht angebbaren Größen ®, N 
und ©: 


Das Integral in (17) läßt sich leicht ausführen und die dritte Näherung bekommt, wenn man 
F(3) nach dem beim Kreisbogen erörterten Konstruktionsprinzip ermittelt, dann die Gestalt: 


3 


3 45° a? 3° (3° — 


a? 
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Führen wir jetzt den Parameter er 4 ein, so erhalten wir die dritte Näherung, 
indem wir bei der Entwicklung nach e bis zum Gliede in gehen. Die Formel (33) zeigt 
nun aber darüber hinaus, daß y” — u K“® für e=0 nicht mehr analytisch in e ist. Vielmehr 
folgt auf das Glied mit e’ ein Glied in e*Ine usw.°): 

3) 


Auf den dritten Schritt hat dies an sich noch keinen Einfluß, da für die Herleitung der 
Differentialgleichungen (2a), (2b), (2c) bzw. (3a), (3b), (3c) nur dreimalige Differenzierbarkeit 
von y nach A gebraucht wurde. Für die weiteren Schritte ist dann natürlich aber eine be- 
sondere Sorgfalt nötig. Hier sei nur noch darauf hingewiesen, daß w“’, wenn man den An- 
satz (1) für y dahin erweitert, daß man auf das Glied in e* einen Term in e*Ine folgen läßt, 
die Vorzahl von e'Ine | angibt. Der Abschnitt 


von y genügt nämlich der Gl. (3c). Genau derselben Differentialgleichung genügt aber auch 
der Abschnitt: 


Die Vorzahl k,(z,y) von 4*InA erfüllt die Gleichung: 
+k,yy=0. 
Da wir (3c) durch (4c) ersetzen und für y“® die genauen Benindingingen von y vor- 


schreiben, wird also durch y“” auch k, (x, y) richtig erfaßt. 


Wir werden also bei der Auswertung in der Entwicklung nach e& bis eine gehen. 
(Beim Beispiel A des affın verzerrten Kreisbogens hatten wir uns schon bei der Konstruktion 
von K“” auf die Glieder bis f? beschränkt. Der geschlossene Ausdruck von K“” bleibt ana- 


lytisch in f.) 

Bild 8 zeigt die Maximalgeschwindigkeit der Ellipse der Dicke 5 wie sie sich mit der 

aus (33) berechneten Formel bis zu den Gliedern in «*Ine ergibt: 

(5 

3 


| 


7 7 3 


x 


9 x 1 A. 1 /x? 
| 


®) Entwickelt man die Rayleighsche Näherung nach einem Parameter 4, so kommen u. U. solche logarith- 
mische Glieder vor. Daß sie auch bei einem an die Prandtlsche Näheruug anschließenden Verfahren auftreten, 


wurde von C. Schmieden bemerkt. 


wobei Ge 
a® 3+ a? 
Die 
Ra 
(34). 

= Bil 
+e*(Ine} 
wi 
er\ 
un 
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Gegenüber der zweiten Näherung hat man eine weitere Erhöhung der Übergeschwindigkeit. 
Zum Vergleich ist noch die Rayleighsche’) Näherung eingetragen. Beschränkt man sich 
in (34) auf die Glieder bis zum Quadrat der Machschen Zahl a, so erhält man: 


Dieser Ausdruck stimmt mit der nach e entwickelten Maximalgeschwindigkeit nach der 
Rayleighschen Näherung überein. 


93 
[7 


Bild 8 (links). Maximalgeschwindigkeit bei der 
nichtangestellten Ellipse der Dicke in Abhängig- 


keit der Machschen Zahl. 
- Kurve.l, 2, 3: erste, zweite und dritte Näherung. 
aa 4: Rayleigh 2. Näherung (Kaplan). 
Anströmgeschwindigkeit 
a. Schallgeschwindigkeit } Im Unendlichen. 
Umax Maximalgeschwindigkeit. 


N pe: 


45 05 5 


Bild 9. Die kritische Anströmgeschwindigkt it Z Bild 10. Die kritische Anströmgeschwindigkeit z 
als Funktion des Achsenverhältnisses der Ellipse. als Funktion des Achsenverhältnisses der Ellipse. 


Kurve 1, 2, 3: erste, zweite und dritte Näherung. Kurve 1: inkompressibel. 

U Anströmgeschwindigkeit Kurve 2: Rayleigh 2. Näherung (Kaplan). 

a% Schallgeschwindigkeit } im Unendlichen. Kurve 3: Tsien. 

a* kritische Schallgeschwindigkeit. Kurve 4: dritte Näherung. 
U Anströmgeschwindigkeit 
a, Sehallgeschwindigkeit } im Unendlichen. 
a* kritische Schallgeschwindigkeit. 


Der Eintritt der kritischen Geschwindigkeit wird in Bild 8 durch die Kurve € angegeben. 
Bild 9 zeigt die kritische Anströmgeschwindigkeit als Funktion des Dickenverhältnisses. 
Wenn auch der ganze Kurvenverlauf bis e=1 gute Werte annimmt, waren, da wir eine Ent- 
wicklung nach der Dicke vorgenommen haben, an sich nur für kleine e. richtige Werte zu 
erwarten. 


Bild 10 gibt einen Vergleich unserer Näherung mit den von Kaplan’) nach Rayleiglı 
und Tsien®) nach seiner Methode erhaltenen Werten. 484 


?,C. Kaplan: Two dimensional subsonie compressible flow past elliptie eylinders, NACA Report 624. 
s) H. S. Tsien: Two-dimensional subsonie flow of compressible fluids. J. Aeron. Sei. Bd. 6 (1939, S. 399 bis 407. 
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Der durch Kapillarwirkung bedingte Zusammenhalt zweier 
benetzter Körper. 
Von Friedrich Schiel in Dresden. 
Mitteilung aus dem Erdbauinstitut der Technischen Hochschule Dresden. 


Ausgehend von der bekannten Haftungskraft zweier benetzter Kugeln gleicher 
Größe wird der Zusammenhalt beliebig geformter Körper untersucht. Für 
den Grenzwert bei verschwindender Benetzungsflüssigkeit ergeben sich be- 
merkenswerte Zusammenhänge mit einigen Begriffen der Differential- 
geometrie. 


1. Einleitung. 


Der Zusammenhalt feuchten Sandbodens entsteht durch die Kapillarkräfte des an den 
Berührungsstellen der Körner befindlichen Wassers. Zur Erklärung dieser sogen. „schein- 
baren Kohäsion“ berechnete Versluys') die Kraft, mit der zwei gleichgroße Kugeln zu- 
sammengehalten werden, wenn an ihrer Berührungsstelle sich Wasser befindet, das die Kugeln 
vollständig benetzt (Randwinkel Null). Als kapillare Oberfläche nahm er eine Drehfläche mit 
kreisförmigem Meridian an, die sich berührend an die Kugeln legt. Er fand, daß die Haftungs. 
kraft dem Kugelhalbmesser R verhältnisgleich ist und daß diese Kraft beim Verschwinden 
des Wassers durch Austrocknen dem Grenzwert 2x0 R zustrebt, worin o die Oberflächen- 
spannung des Wassers ist. In einer Kugelschüttung ist die Anzahl der Kugeln auf der Flächen- 
einheit 1/ R? verhältnisgleich, so daß die auf die Flächeneinheit bezogene „scheinbare Kohäsion“ 
einer feuchten Kugelschüttung mit 1/R verhältnisgleich zunimmt, was mit dem Verhalten 
feuchten Sandbodens durchaus im Einklang steht. 


Im Anschluß an die Ergebnisse von Versluys erheben sich einige Fragen, die im 
folgenden behandelt werden sollen. 


1. Ist die Annahme kreisförmigen Meridians für die kapillare Drehfläche zulässig oder 
ergibt sich bei Einführung der richtigen, nach der Kapillaritätstheorie zu bestimmenden Form, 
ein anderes Bild? 


2. Strebt bei unvollständiger Benetzung die Haftungskraft ebenfalls einem endlichen 
Grenzwert zu, oder ist er Null? Man beobachtet nämlich beim Austrocknen feuchter Sand- 
körper einen Zerfall, ehe die letzten Wasserspuren verschwunden sind, was mit einem Null- 
werden des Grenzwertes gut übereinstimmen würde. 


3. Wie groß ist die Haftungskraft bei zwei beliebig geformten, benetzten Körpern, die 
sich in einem Punkt berühren ? 


4. Werden außer der Haftungskraft auch andere (drehende) Kräfte durch das Wasser 
erzeugt? Das Vorhandensein derartiger Kräfte wird durch die Beobachtung nahegelegt, daß 
beim Benetzungsvorgang an der jeweiligen kapillaren Oberfläche ruckartige Bewegungen ein- 
zelner Körner stattfinden, die beim Auf- und Absteigen des Wassers jedesmal eine Setzung 
der Sandmasse bewirken, welche erst nach oftmaliger Wiederholung der Benetzung abklingt?). 


2. Kapillare Drehflächen bei Vernachlässigung der Schwere. 


Der Einfluß der Schwerkraft auf die Veränderung des Flüssigkeitsdruckes p mit der 
Höhe kann bei den hier beabsichtigten Untersuchungen vernachlässig’ verden. Die Kapillaritäts- 
theorie lehrt, daß der Flüssigkeitsdruck nach dem ersten Satz von Laplace der mittleren 
Krümmung der Flüssigkeitsoberfläche verhältnisgleich ist. Bei konstantem p hat so- 
mit die kapillare Oberfläche konstante mittlere Krümmung. Die Drehflächen dieser Eigen- 
schaft sind wohl bekannt und genau untersucht worden’). 

Wir betrachten hier die kapillaren Drehflächen im Hinblick auf die Kraft, welche die 
von ihnen begrenzte Flüssigkeit zwischen zwei festen Körpern übertragen kann Aus Gleich- 
gewichtsgründen muß sich in jedem Schnitt durch die Flüssigkeit durch Zusammensetzung 
der im Schnitt wirkenden Flüssigkeitsdrücke und Oberflächenspannungen dieselbe Kraft er- 
geben. Diese Bedingung führt ebenso zur Bestimmung der Form kapillarer Oberflächen wie 
die Forderung konstanter mittlerer Krümmung; die beiden Bedingungen sind identisch. 


1) De capillaire werkingen in den bodem. Amsterdam 1916, 


R 2) Über derartige Benetzungsversuche berichten Kahl, Manz und Neumann in: Bautechn. Bd. 19 (1941), S. 349 
bis 352, x 


3) Siehe z. B. Plateau: Statiqne experimentale et theorique des liqynides, chap. 11., Gand, 1873. 
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Es werden hauptsächlich folgende Bezeichnungen verwendet: 


p — Flüssigkeitsdruck. 
p;= — p = Zugspannung der Flüssigkeit. 
Oberflächenspannung der Flüssigkeit. 
rt; — Hauptkrümmungshalbmesser der Flüssigkeitsoberfläche, positiv nach dem Flüssig- 
keitsinneren. 
r — Halbmesser einer Flüssigkeitskugel. 
pP = Kraft, die eine ‘von einer kapillaren Drehfläche begrenzte Flüssigkeitsmasse in 
Richtung ihrer Achse überträgt. i 
x|z — Ebenes Koordinatensystem zur Darstellung des Meridians der kapillaren Dreh- 
fläche (z in Richtung der Achse). 
= & in — Halbmesser des „Halsquerschnittes“ der Drehfläche. 
x 
= Krümmungshalbmesser des Meridians an der Stelle x, (Halsquerschnitt). 
Weitere Bezeichnungen werden an geeigneter Stelle erklärt. 
Nach Laplace ist 
p=o\,t7,. (1) 
a) Kugel. Der Druck in einer Flüssigkeitskugel ist nach Gl. (1): 
20 


” In einem Schnitt nach Bild 1 wird folgende Kraft übertragen: 


Die Kugel ist die einzige kapillare Drehfläche, die ohne Anlehnung an 
feste Körper bestehen kann, da die von einer Kugel begrenzte Flüssig- 
Bild 1. Flüssigkeitskugel. Keitsmasse keine Kraft nach außen überträgt. 


b) Katenoid. Die Drehfläche mit der mittleren Krümmung Null ist das sogenannte 
Katenoid mit der Kettenlinie als Meridian (Bild 2). In einer katenoidförmigen Flüssigkeits- 
masse herrscht der Flüssigkeitsdruck p=0 wie unter einer ebenen Oberfläche. Die Kraft, 
die durch ein Katenoid in seiner Achsenrichtung übertragen wird, errechnet sich nur aus o 
am Umfang des Halsquerschnittes, da p=0 ist: 


Jeder Schnitt senkrecht zur z- Achse liefert 
denselben Wert von P, wie aus den Eigen- 
schaften der Kettenlinie leicht bewiesen 
werden kann. 

Taucht man zwei sich berührende 
feste Kugeln in eine benetzende Flüssigkeit 
und hebt sie wieder heraus, so verbleibt 
um die Berührungsstelle ein Flüssigkeits- 
ring nach Bild 3 bestehen, der um so ge- 
nauer katenoidförmig ist, je kleiner die 


Kugeln sind, je weniger also die Schwer- Bild 2. Katenoid. 
kraft ausmacht. Aus der Bedingung, daß 
die Kettenlinie als Meridian die Kugeln Bild 3 (rechts) Zwei Kugeln 


berühren muß, läßt sich errechnen: mit Katenoid. 


x, = 0,6815 R | 
P=1363r0R 
Setzt Verdunstung ein, so wandert die kapillare Oberfläche nach innen, p wird negativ und 
P nimmt zu, wie sich aus einer nach Versluys durchgeführten Näherungsrechnung mit 
kreisförmigem Meridian ergibt. 
c) Nodoid. Die Gleichung des Meridians für Drehflächen mit Flüssigkeitszug p, = —p 
lautet ®): 


(3). 


4) Siehe z. B. Bakker: Kapillarität und Oberflächenspannung, S. 123. Leipzig 1928 (in: Handbuch der Experi- 
ınentalphysik, hrsg. von Wien und Harms, Bd. 6). Die Gleichung ist erstmalig von Beer 1857 aufgestellt worden. 
14 
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Darin sind F(p,k) und E(,k) die elliptischen Integrale erster und zweiter Gattung mit 
dem Modul 
— 

k 

während x von @ folgendermaßen abhängt: 

Die Form der durch Gl. (4) gegebenen Kurve ist in Bild 4 
zu sehen. Als Meridian unserer Drehfläche kommt das dick 
ausgezogene Stück in Frage. Die gesamte Drehfläche führt 
nach Plateau den Namen Nodoid. 

Die Abszisse des Kurvenpunktes mit einer zur Dreh- 
achse senkrechten Tangente bezeichnen wir mit X. Sie ist 
das geometrische Mittel aus x, und «,. Führt man den 
Modul in der üblichen Art ein: 


| ame k=sina, 
Bild 4. Meridian des Nodoids e R 
für a—80°, so ergibt sich aus den gegebenen Formeln 
x x 
—=cosa —=ycosa. 
2 x 


Die Kraft, die durch das Nodoid übertragen wird, berechnen wir einmal durch die 
Kräfte im Halsquerschnitt 


und durch einen Schnitt an der Stelle e=x, an der o keinen Beitrag liefert, 
P=n2p»,. 


Durch Ausschaltung von p, und Einführung der obigen Beziehungen erhält man 


Je mehr sich der Modul k dem Wert 1, bzw. a dem Wert 5 nähert, desto mehr nähert sich 
der Meridianbogen der Gestalt einer Kettenlinie, und für «a => geht (5) in Gl. (2) über, die 
für P beim Katenoid gilt. 

Um zu sehen, welche Form der Meridian des Nodoids bei klein werdendem a annimmt, 
berechnen wir den Krümmungshalbmesser o an der Stelle <=x,. Nach Gl. (1) ist (o positiv 
einzuführen) 

2x, 


mit Rücksicht auf die beiden obigen Gleichungen für P. Daraus erhält man 


%, 


Beim Grenzübergang lim a—=0 nähert sich 


dem Wert 1 und o dem Wert 
2 


3 Mit anderen Worten: der Meridian 


nähert sich der Kreisform. Wegen der 
Annäherung von x, an x, bei lima=0 
wird dabei o klein gegen x,. In Bild 5 
ist der Meridian für a=45° gezeichnet, 
und man erkennt, daß er von der Kreis- 
form selbst bei diesem großen Winkel nur 
noch wenig abweicht. 

Die Kraft P wird für a=0 im allgemeinen unendlich groß. Wenn gleichzeitig x,=0 
wird, ist P unbestimmt. Welchen Grenzwert P annimmt, hängt davon ab, wie die beiden 


Bild 5. Nodoid mit @«—=45®, 
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Grenzübergänge lima=0 und limx,=0 miteinander verbunden sind, und man kann nach 
obigem zur Berechnung des Grenzwertes von P die Meridianlinie kreisförmig an- 
nehmen, sofern es sich um einen Meridian in einem engen Spalt handelt (o klein gegen «,). 

Die Drehflächen mit positivem p (Unduloid und äußerer Teil des Nodoids) werden hier 
nicht benötigt und daher übergangen. 


3. Grenzwert der Haftungskraft zweier benetzter Kugeln. 


a) Vollständige Benetzung. In der eingangs erwähnten Arbeit hat Versluys für die 
Haftungskraft P zweier Kugeln eine Formel aufgestellt, die für das Verschwinden der be- 
netzenden Flüssigkeit einen unbestimmten Ausdruck liefert, der sich der Grenze 2r Ro nähert. 
Man kann diesen Grenzwert unmittelbar erhalten (siehe Bild 6). Wir bezeichnen den Halb- 
messer des Halsquerschnittes wieder mit «, und mit 2 die Ordinate der Benetzungsgrenze 
nach Bild 6. Wenn x, klein gegen R und somit z klein gegen «,- ist, gilt näherungsweise 


P=nap, +t2n2,0=2noR+2nox,, 
Bild 6. Zwei Kugeln. 


Durch die vollständige Verdunstung wächst die Zugspannung der Flüssigkeit von 0 
auf ©, und die Haftungskraft nimmt von dem Wert nach Gl. (3) auf jenen nach Gl. (6), also 


“um knapp zu. 


b) Unvollständige Benetzung. Bildet die Flüssigkeit mit den festen Kugeln einen Rand- 
winkel #, so läßt sich der Grenzwert von P durch eine ähnliche Näherungsrechnung erhalten 
(siehe Bild 7) 

220:.c08sd 


Die Näherung ist zulässig, solange # < 3 und wir beschränken 
uns auf dieses Intervall. 


_0 _20oRcosd Bild 7. Unvollständige Benetzung. 
2 
x 
x 


Der Grenzwert ist nicht Null, sondern endlich, solange 9 < 3 ist. Das Zerfallen von feuchten 


Sandkörpern vor restloser Austrocknung kann also nicht an unvollständiger Benetzung liegen, 
sondern muß andere Ursachen haben. 


4. Beliebiger Körper und Berührungsebene. 


a) Grenzwert der Haftungskraft. Bei der Ableitung der Gl. (6) für P bei zwei Kugeln 
konnten wir bemerken, daß infolge des Unendlichwerdens von p, das Glied mit o neben 
jenem mit p, vernachlässigt werden konnte. Aufbauend auf dieser Feststellung läßt sich 
eine entsprechende Näherungsrechnung für die Haftungskraft anstellen, welche beim Ver- 
dunsten einer Flüssigkeitsmenge entsteht, die den Berührungspunkt eines beliebigen Körpers 
mit seiner Berührungsebene umgibt. 

Wir setzen voraus, daß die Oberfläche des Körpers an der Berührungsstelle stetig und 
bauchig gekrümmt oder der Berührungspunkt ein elliptischer Flächenpunkt ist, schalten also 
parabolische Flächenpunkte aus, während sich sattelförmige Krümmung oder hyperbolische 
Flächenpunkte von selbst erübrigen. 

Die Abstände der Benetzungsgrenze von der Berührungsebene bezeichnen wir wieder 
mit z und legen in diese Ebene ein Koordinatensystem, dessen Achsen x/y mit den Haupt- 
krümmungsrichtungen der Körperoberfläche zusammenfallen. Außerdem verwenden wir 
folgende Bezeichnungen: 


Nx,N0, Krümmungen (1/r) der Hauptnormalschnitte der Körperoberfläche, 


&%,,Y, =kleinste Entfernungen der kapillaren Oberfläche von der Flächennormalen, ge- 
messen in Richtung der Achsen. 


14® 
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Bei fortschreitender Verdunstung des Flüssigkeitsrestes 
nähert sich die Projektion des Halsquerschnittes auf die 
Berührungsebene der Gestalt der Dupinschen Kenn- 
figur°). Der Abstand z (Bild 8) ist beim Grenzübergang 
längs des ganzen Umfanges der Kapillarfläche konstant und 
es ist näherungsweise 


Bezeichnen wir die Fläche der Ellipse mit den Halbachsen 

x%,, y, mit J, so ist demnach 

yK Dupinschen Kennfigur. 

worin K=n,.n, das Gaußsche Krümmungsmaß der Körperoberfläche ist. Auf dieser Fläche J 

wirkt die Zugspannung 


J 


und man erhält 4no 


YK . . . . . . . . . . . . . (9). 


Es ließe sich hier ebenso wie bei den zwei Kugeln eine Näherungsrechnung für den 
Fall unvollständiger Benetzung anstellen, doch wollen wir uns hier und im folgenden aus 
Raumgründen auf den wichtigsten Fall der vollständigen Benetzung beschränken. 


b) Gleichgewichtslagen des Körpers. Der Körper kann sich in bezug auf die Kraft P nur 
dann im Gleichgewicht befinden, wenn P in Richtung der Flächennormalen durch den Be- 
rührungspunkt geht. Dies wird um so genaucr zutreffen, je kleiner die Flüssigkeitsmenge ist, 
d. h. je mehr sich ihre Form der Dupinschen Kennfigur annähert. Bei endlicher Größe der 
Flüssigkeitsmenge liegt ihre äußere Begrenzungsfläche nicht genau symmetrisch zur x- und 
y-Achse, so daß die Kraft P nicht genau im Berührungspunkt angreift. Infolgedessen hat 
der Körper das Bestreben, auf seiner Berührungsebene abzurollen. Das Drelımoment, das den 
Körper dabei antreibt, wird offenbar beim Verdunsten mit z2 zusammen gleich Null, strebt 
also nieht einem endlichen Grenzwert zu wie die Haftungskraft. 

Das größtmögliche abrollende Moment läßt sich nicht in allgemeingültiger Form an- 
geben, da es von der gesamten Oberflächenform des Körpers abhängt. Wir wollen daher die 
Frage nach dem abrollenden Moment in folgender beschränkter Form lösen: 1. in welcher 
Richtung rollt der Körper und 2. in welcher Lage befindet er sich im Gleichgewicht? Wir 
nehmen an, daß die Flüssigkeitsmenge V, die den Berührungspunkt umgibt, beim Abrollen 
unverändert bleibt. Es läßt sich zeigen, daß bei kleinem 2 die mittlere Entfernung der 
Körperoberfläche von ihrer Berührungsebene für das ganze Gebiet der Kennfigur gleich der 
halben Randentfernung 2 ist. Dementsprechend ergibt sich die Flüssigkeitsmenge V 


Schaltet man aus dieser Gleichung und aus (8) den Abstand z aus, so erhält man: 

Der Körper wird so abrollen, daß sich die freie Oberflächenenergie «er Flüssigkeit ver- 
mindert, und wird dann eine Gleichgewichtslage einnehmen, wenn diese Energie einen 
extremen Wert hat. Bezeichnen wir mit F die gesamten, der Flüssigkeit zugänglichen festen 
Oberflächen (abrollender Körper und Berührungsebene), so ist die gesamte Flüssigkeits- 
oberfläche 

F—-2J, 

da sowohl auf dem Körper als auch auf der Berührungsebene von der gesamten benetzten 
festen Oberfläche ein Stück von der Größe J fehlt. Dabei ist das ringförınige Stück Flüssig- 
keitsoberfläche von der Höhe z längs des Umfanges der Kennfigur gegenüber J vernachlässigt, 
was bei der vorausgesetzten Kleinheit von z zulässig ist. 

Die freie Oberflächenenergie ist somit 

E=os(F—-2J). 


5, Diesen Verdeutschungsvorschlag für „Indikatrix“, sowie wertvolle Ratschläge hinsichtlich der im folgenden 
verwendeten Benennungen neuer Begriffe verdanke ich meinem Bruder Theodor Schiel in Berlin, Durch die Zu- 
sammenstellung des Wortes „Kennfigur‘* mit „Dnpin“, dem Namen des Urhebers, ist wohl ohne weitere Erklärung ver- 
ständlich, was gemeint ist. 
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Beim Abrollen um ein Stück ds ist das Energiegefälle 


Aus Gl. (10) folgt: 
dJ _YaV dK 
ds 2yK® ds’ 
dK 
de de 


Den Ausdruck ee kann man das Krümmungsgefälle nennen. Der Körper rollt in 


jener Richtung ab, inder dasKrümmungsgefälle am größten ist. Das Energie- 


gefälle liefert dimensionsmäßig nicht ein Moment, sondern eine Kraft, so daß statt des oben 
gebrauchten Ausdruckes „abrollendes Moment“ genauer „abrollende Kraft“ zu sagen wäre. 

Um sich ein vollständiges Bild vom Verhalten des Körpers zu machen, muß man auf 
seiner Oberfläche jene Punkte verbinden, in denen das Gaußsche Krümmungsmaß_ den- 
selben Wert hat. Es ergibt sich die Schar der „Linien gleichen Krümmungsmaßes“. Senk- 


recht zu diesen ist das Krümmungsgefälle am größten. Die verschiedenen möglichen, jeweils 


in Richtung des größten Krümmungsgefälles verlaufenden Linien bilden eine weitere Schar; 
man kann sie die „Krümmungsfallinien“ nennen. Die Krümmungsfallinien gehen von den 
Krümmungsmaxima aus und vereinigen sich in den Krümmungsminima. 

Ähnlich wie eine schwere Kugel auf einer schiefen Ebene in Richtung der „Fallinie“ 
abrollt, wird ein gewichtslos gedachter benetzter Körper, der mit irgendeinem Punkt seiner 
Oberfläche auf eine benetzte Ebene aufgesetzt wird, längs seiner durch diesen Punkt gehenden 
Krümmungsfallinie abrollen. Er wird erst zur Ruhe kommen, wenn sein Berührungspunkt 
mit einem Krümmungsminimum seiner Oberfläche zusammenfällt, wobei er eine stabile Gleich- 
gewichtslage erreicht hat. Setzt man ihn mit einem Krümmungsmaximum auf, so befindet 
er sich im labilen Gleichgewicht. Man kann also kurz sagen, die benetzende Flüssigkeit ver- 
anlaßt den Körper, sich „möglichst flach“ auf die Ebene zu legen, was sich am Verhalten 
kleiner wasserbenetzter Körper leicht bestätigen läßt. Bei größeren Körpern müßte man den 
Einfluß der Schwerkraft z. B. dadurch ausschalten, daß man außer der benetzenden Flüssig- 
keitsschicht eine zweite Flüssigkeit verwendet, die den ganzen Körper umgibt und seine 
Schwere durch Auftrieb aufhebt. 

Bei Drehflächen sind die Parallelkreise Linien gleichen Krümmungsmaßes und die Me- 
ridiane Krümmungsfallinien. Bei Drehflächen sind also die Krümmungsfallinien gleichzeitig 
Krümmungslinien, was bei beliebigen Flächen nicht der Fall ist. Drehflächen haben außer- 
dem die Eigentümlichkeit, daß die Krümmungsminima 
meist längs eines ganzen Parallelkreises auftreten. Ein 
spindelförmiger Drehkörper nach Bild 9, oben, ist daher u 
abgesehen von der Schwerkraft immer im indifferenten + 
Gleichgewicht, wenn er mit seinem größten Parallelkreis 
auf eine benetzte Ebene gesetzt wird. Eine flachgedrückte 
Spindel nach Bild 9, unten, ist dagegen auf der flachen —--2 
Stelle aufsitzend immer im stabilen Gleichgewicht. —— 

Das Bestreben des Körpers abzurollen wächst mit 


der Größe des Krümmungsgefälles („Feldstärke‘). 


5. Zwei beliebige Körper. 


a) Darstellung der Krümmungsverhältnisse am Berührungspunkt. Wir setzen stetige 
Krümmung der beiden Körperoberflächen voraus, wollen uns aber nicht auf elliptische Punkte 
beschränken. Wir setzen nur voraus, daß die beiden Oberflächen sich berühren, ohne sich 
zu durchdringen, was auch bei gegenseitigen Verdrehungen um die gemeinsame Flächen- 
normale zutreffen soll. Auf das Abrollen der beiden Körper aufeinander gehen wir nicht ein. 
Es soll z. B. dadurch verhindert sein, daß sich zwei Krümmungsminima gegenüberliegen. 

Um die Verhältnisse leichter übersehen zu können, bedienen wir uns einer zeichnerischen 
Darstellung des Eulerschen Satzes über die Krümmung der Normalschnitte einer Fläche: 

Zwischen diesen Krümmungen und den Normalspannungen eines ebenen Spannungszustandes 
besteht eine vollständige Analogie (Krümmungstensor — Spannungstensor), so daß man eine 
dem Mohrschen Spannungskreis entsprechende Darstellung der Krümmungen benützen kann. 
Wir nennen diese Darstellung das Krümmungsbild (Krümmungskreis würde zu Mik- 
verständnissen führen). 
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Trägt man (siehe Bild 10) auf einer Achse die beiden Hauptkrümmungen », und n, auf 
und zeichnet über n, —n, als Durchmesser einen Kreis, so läßt sich für jeden Winkel 
eines Normalschnittes mit dem Hauptnormalschnitt %», nach Bild 10 die Krümmung n 
ablesen. Der Abstand des Kreismittelpunktes vom Ursprung der n-Achse ist die mittlere 
n,+n, 


Krümmung m Das geometrische Mittel aus n, und n,, das gleich der Wurzel 


aus dem Gaußschen Krümmungsmaß ist, läßt sich als Länge der Kreistangente ablesen. 
Bild 10a ist die Darstellung eines bauchig gekrümmten (elliptischen) Punktes. Bei para- 
bolischen Punkten (Bild 10b) geht der Kreis durch den Anfangspunkt der Krümmungsachse. 
Bei sattelförmig gekrümmten (hyperbolischen) Punkten läßt sich Y— K nach Bild 10c ablesen. 


n ] positive Rıchtung 
4 Bild 11 (rechts). Berührung der Flächennormalen 


Bild 10. Darstellung des Eulerschen Satzes Bild 12. Krümmungsbild zweier sich berührender 
durch das „Krümmungsbild‘. Flächen. 

Auf der gemeinsamen Flächennormalen der beiden sich berührenden Körperoberflächen 
bezeichnen wir eine Richtung als positiv. Alle Krümmungen, die vom positiven Ende der 
Flächennormalen aus gesehen erhaben erscheinen, seien positiv, die hohlen negativ. Die Krüm- 
mungen der Normalschnitte der einen Fläche bezeichnen wir mit n, die der anderen mit n, 
den gegenseitigen Abstand der Flächen, gemessen senkrecht zur gemeinsamen Berührungs- 
ebene, mit 2. Für einen Punkt im sehr kleinen Abstand x vom Berührungspunkt gilt (siehe 
Bild 11) näherungsweise 


2 
nz. 


Wir definieren 


als die „gegenseitige Krümmung“ der beiden Flächen. 

Für die gegenseitige Krümmnng N, gebildet für alle Richtungen in der Berührungsebene, 
gilt der Eulersche Satz. Es gibt also zwei Extremwerte N,, N, und für die Richtung, die 
mit der Richtung von N, den Winkel @ einschließt, ist 

N=N, cos®p+ N, sin’ (11). 
Es erübrigt sich die Führung eines unmittelbaren Beweises aus dem Bestehen des Euler- 
schen Satzes für n und n, da die Voraussetzungen seiner Herleitung, nämlich die Darstellungs- 
möglichkeit der Abstände einer Fläche von ihrer Berührungsebene durch eine Potenzreihe, 
auch für die gegenseitigen Abstände der beiden Flächen erfüllt ist, wenn voraussetzungs- 
gemäß am Berührungspunkt beide Flächen stetig sind. 

Aus (11) folgt, daß der geometrische Ort gleicher Abstände 2 der beiden Flächen mit 
klein werdendem z die Gestalt eines Kegelschnittes annimmt. Wir nennen ihn die Be- 
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rührungskennfigur der beiden Flächen und unterscheiden dementsprechend ellip- 
tische, parabolische und hyperbolische Berührung zweier Flächen. Bei Berührung ohne Durch- 
dringung ist die Kennfigur eine Ellipse. Aus den beiden gegenseitigen Hauptkrümmungen 
N,, N, berechnen wir das „Maß der gegenseitigen Krümmung“ 

Die Dupinsche Kennfigur und das Gaußsche Krümmungsmaß ergeben sich als Sonderfälle, 
wenn eine der beiden Flächen eben ist. 

Die Berührungskennfigur und die Werte N,, N, und K ändern sich bei gegenseitigen 
Verdrehungen der beiden Körper um die gemeinsame Flächennormale im Berührungspunkt. 
Der Verlauf dieser Größen bei Verdrehungen läßt sich am Krümmungsbild verhältnismäßig 
einfach verfolgen. In Bild 12a wurden die Hauptkrümmungen beider Körper am Berührungs- 
punkt auf dieselbe n- Achse aufgetragen und durch Zeiehnuug der beiden Kreise mit den 


Nn,—N, 


Halbmessern a = und =". das Krümmungsbild erhalten. Die mittleren Krüm- 


mungen der beiden Flächen seien m und m, und wir bezeichnen 
M=m-—m 

als die „mittlere gegenseitige Krümmung“. Die Hauptkrümmungsrichtung #, sei gegenüber n, 
um den Winkel » verdreht, und der Normalschnitt x — x (siehe Bild 12b), für den wir die 
gegenseitige Krümmung N bestimmen wollen, schließe mit n, den Winkel 9, mit ”, den 
Winkel 9+® ein. Die Krümmungen n und » können für diese Schnittrichtung aus Bild 12a 
entnommen werden, und es ist N=n—n. Man erkennt aus Bild 12c, daß N gleich der 
Summe von M und der Projektion der Strecke AB auf die n-Achse ist. Läßt man bei fest- 
gehaltenem w (feste Lage der beiden Körper) p wachsen (Drehung der Normalschnittebene), 
so durchläuft diese Projektion bei Drehung des Dreieckes A BC den Bereich — AB bis + AB. 


Demnach ist Na=N=M+AB, 
Nyn=N,=M—-AB 
und das Maß der gegenseitigen Krümmung ist 
K=NN,=M'’— AB, 
Wenn sich die beiden Hauptnormalschnittebenen beider Körper decken, also für o=0 und 


= 3 nimmt K extreme Werte an. Die Darstellung dieser Werte im Krümmungsbild wird 
im nächsten Abschnitt erörtert. 


b) Grenzwert der Kräfte im Berührungspunkt. Bei Betrachtung der kapillaren Dreh- 
fläche konnten wir davon ausgehen, daß die in jedem Schnitt durch diese Fläche an- 
greifenden Flüssigkeitsdrücke p und Oberflächenspannungen o als Resultierende stets dieselbe 
Kraft P ergeben müssen. Im vorliegenden allgemeinen Fall muß an die Stelle der Kraft P 
ein räumlichesKraftkreuz treten, das wir durch eine Kraft P in Richtung der Flächen- 
normalen und ein um diese Linie drehendes Drehmoment D bestimmen. 

An der Berührungskennfigur lassen sich zur Berechnung von P dieselben Überlegungen 
anstellen wie an der Dupinschen Kennfigur bei Herleitung der Gl. (9). Es gilt daher un- 


verändert 


nur bedeutet hier K nicht das Gaußsche Krümmungsmaß, sondern das Maß der gegenseitigen 
Krümmung. 

Eine genaue Ermittlung des Drehmomentes D bei gegebener Form der beiden sich be- 
rührenden Körper ist hier nicht beabsichtigt, sondern es sollen die Werte von D ermittelt 
werden, die sich bei fortschreitender Verdunstung der Flüssigkeit einstellen. Wir machen 
ganz entsprechende Näherungsannahmen wie bei der Behandlung des Abrollens eines Körpers 
auf einer Ebene: Die Flüssigkeitsmenge sei so gering, daß sich ihre Form der „Berührungs- 
kennfigur“ genügend annähert, und daß die freie Flüssigkeitsoberfläche, die sich in Höhe des 
gegenseitigen Flächenabstandes z ausbildet, gegenüber der Fläche der Kennfigur vernachlässigt 
werden kann. Weiter nehmen wir an, daß bei gegenseitigen Verdrehungen der beiden Körper 
die Flüssigkeitsmenge V erhalten bleibt. 

Der Flächeninhalt J der Berührungskennfigur läßt sich durch ganz entsprechende Über- 
legungen finden wie bei der Herleitung von Gl. (8), und es ergibt sich 
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worin K das Maß der gegenseitigen Krümmung ist. Für die Flüssigkeitsmenge gilt wie bei 
der Dupinschen Kennfigur 


| 
Wir berechnen wieder die gesamte freie OÖberflächenenergie: 
(16). 


Bei einer gegenseitigen Verdrehung beider Körver um den Winkel & ändert sich E, und das 
durch die Flüssigkeit ausgeübte Drehmoment ist gleich dem Energiegefälle: 


dE dJ dJ_ aVdK 


Das Drehmoment D ändert sich periodisch mit & und ist an den Stellen ®=0 und 
o=5 gleich Null, an denen K und damit P nach Gl. (14) Extremwerte hat. Für diese 


Stellen können die Werte von K, die der 
Haftungskraft P umgekehrt proportional sind, 
nach Bild 13 als geometrische Mittelwerte 
von N, und N, aus dem Krümmungsbild ent- 


nommen werden. Die stabile Gleichgewichts- 
durch gek ichnet, daß di 
lage ist dadurch gekennzeichnet, da ie 
Werte N, und N, die größtmögliche Ver- Bild 13. Krümmungsbild für die Gleichgewichtslagen. 
schiedenheit aufweisen, oder die Exzentrizität Krümmungsmaß: Kyax| 
i Haftungskraft: Pmin | P max» 
der Kennfigur Ellipse am größten ist. Handelt | 
es sich um zwei bauchig gekrümmte Körper, Oberflächenenergie: Emax Enmin» 
die sich mit ihren erhabenen Seiten berühren, Gleichgewicht: labil' stabil. 


so decken sich in der stabilen Gleichgewichtslage die Richtungen geringster Krümmung. 
Ist der Berührungspunkt auf einer der beiden Körperoberflächen Nabelpunkt (n,—n,), so 
ist D beständig Null, wie sich aus (17) ergibt. 

Für den Fall, daß die Stetigkeitsbedingung nicht erfüllt ist und der eine Körper am 
Berührungspunkt eine Spitze hat, läßt sich leicht beweisen, daß die Haftungskraft beim Ver- 
schwinden der Flüssigkeit dem Wert Null zustrebt (K=o in (14)). Diese Tatsache kann 
für das erwähnte vorzeitige Zerfallen feuchter Sandkörper mit verantwortlich sein, da infolge 
der Kristallstruktur selbst scheinbar ganz rundliche Sandkörner feine Ecken und Kanten haben. 


c) Zusammenfassung. Bringt man zwei stetig gekrümmte Körperoberflächen, die mit 
einer benetzenden Flüssigkeitsschicht bedeckt sind, in Berührung, so werden durch die Ober- 
flächenspannung Kräfte hervorgerufen, die zwischen den beiden Körpern wirken. Zu ihrer 
kurzen Kennzeichnung wurde der Begriff der „Berührungskennfigur“ zweier Flächen eingeführt, 
die, entsprechend der an einer berührenden Ebene auftretenden Dupinschen Kennfigur, als 
geometrischer Ort gleicher Flächenabstände definiert ist. 

Die Kapillarkräfte suchen die beiden Körper so aufeinander abzuwälzen und zu drehen, 
daß eine gegebene, den Berührungspunkt umgebende kleine Flüssigkeitsmenge V eine mög- 
lichst große Fläche J in Form der Berührungskenntigur bedeckt — mit anderen Worten: die 
Körper legen sich „möglichst flach“ aneinander. Die Haftungskraft, die mit abnehmender 
Krümmung zunimmt, ist in der stabilen Gleichgewichtslage am größten. Ist der eine Körper 
eben oder kugelig gekrümmt (Nabelpunkt!), so findet kein Drehen, sondern nur Abwälzen 
statt. Es erfolgt längs der „Krümmungsfallinien“ des anderen Körpers in der Richtung der 
stärksten Abnahme des Gaußschen Krümmungsmaßes. 

Wenn die den Berührungspunkt umgebende Flüssigkeitsmenge abnimmt und allmählich 
verschwindet, strebt die Haftungskraft einem endlichen Grenzwert zu, während die abwälzenden 
und drehenden Kräfte Null werden. Alle diese Kräfte teilen mit den übrigen kapillaren Er- 
scheinungen die Eigenschaft, bei abnehmender Körpergröße (= Korngröße im Sand) gegenüber 
Massenkräften an Bedeutung zuzunehmen. 

Die gefundenen Erscheinungen kann man dazu verwenden, um die Krümmungslinien 
einer Fläche wiederzugeben. Man kann dazu einen Körper nach Bild 9, unten, benützen. Legt 
man diesen „Krümmungskompaß®)“ auf eine benetzte Fläche, so dreht er sich in die Richtung 
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des einen Hauptnormalschnittes der Fläche. Denkt man sich die Fläche als sehr dünne 
Schale körperlich verwirklicht und legt auf den gegenüberliegenden Punkt der benetzten Rück- 
seite den Krümmungskompaß, so zeigt er die Richtung des anderen Hauptnormalschnittes an. 
Das „Direktionsmoment“ nimmt mit zunehmender Angleichung der beiden Hauptkrümmungen 
aneinander ab, und an Nabelpunkten bleibt der Krümmungskompaß in jeder Lage in Ruhe. 474 


Zur Triebstockverzahnung. 
Von Walter Wunderlich in Kiel. 
Die bei der Triebstockverzahnung mit auf dem Teilkreis angeordneten 
Triebstockmittelpunkten auftretende Unstimmigkeit wird gestalts- und 
zahlenmäßig untersucht. Der an sich geringe Fehler kann durch eine 
kleine Verschiebung der Triebstöcke gegen die Radmitte behoben werden. 
Die Mindestverschiebung wird berechnet und ebenso wie die Fehlerschranke 
nomographisch dargestellt. 


Im Rahmen des Kapitels „Verzahnungen“ taucht in den meisten Lehr- und Handbüchern 
der Kinematik und Getriebelehre, wie auch in den Vorlesungen immer wieder das Beispiel 
der Triebstockverzahnung auf, trotz ihrer untergeordneten Bedeutung, wohl wegen 
ihrer offensichtlichen Einfachheit. Dabei ist gerade hier eine grundsätzliche Schwierigkeit 
vorhanden, die — bewußt oder unbewußt — fast stets übergangen wird, da die übrigen 
Sonderformen (Evolventen- und Zykloidenverzahnnungen) nichts Derartiges aufweisen. Im vor- 
liegenden Aufsatz wird die bei der üblichen Annahme zwangläufig auftretende Unstimmigkeit 
der Triebstockverzahnung zahlenmäßig untersucht. Gleichzeitig wird gezeigt, daß der an sich 
sehr kleine Fehler durch eine geringe Abänderung gänzlich behoben werden kann. 


$ 1. Der Fehler der Triebstockverzahnung. 

Bei der Triebstockverzahnung sind bekanntlich auf dem Umfang des einen Rades in 
regelmäßigen Abständen achsenparallele zylindrische Zapfen oder Rollen — die Triebstöcke — 
angeordnet, die im Verlauf der Drehung in die entsprechenden Zahnlücken des anderen Rades 
eingreifen. Es sind demnach die Zahnflanken für das erste Rad in Kreisform p, vorgegeben, 
wodurch das Zahnprofil p, des zweiten Rades bestimmt ist und nach einem der üblichen 
allgemeinen Verfahren konstruiert werden kann. Man kann etwa den ersten Wälzkreis k, 
auf dem zweiten abrollen lassen und hierbei die Hüllkurven p, der mitgenommenen Kreise p, 
ermitteln, wobei es natürlich nur auf die inneren Äste ankommt. 

Da die Mittelpunkte M der Triebstöcke fast ausnahmslos auf dem Wälzkreis ange- 
nommen werden, beschreiben sie bei der Rollung Zykloiden m, und die Gegenprofile p, 
sind deren Parallelkurven im Abstand s, wenn s den Halbmesser von p, bezeichnet 
(Bild 1). In dem Augenblick, in dem der Mittelpunkt M auf dem festen Teilkreis k, aufsetzt, 
befindet er sich in einer Spitze seiner Bahn, die Zahnflanke p, nimmt im Punkt P auf der 
Berührungstangente von k, und k, ihren Anfang, während die symmetrische Flanke p, im 
Gegenpunkt P zu Ende ist. Zwischen P und P fügt sich als Verbindung zwanglos ein 
Halbkreis vom Halbmesser s ein, der die Zahnlücke schließt. Da er sich in P und P 
berührend anschließt, scheint alles in bester Ordnung zu sein '). 

Und doch ist dies nicht 
der Fall, wenn man die Stelle P 
etwas genauer unter die Lupe 
nimmt. Das Profil p, hat näm- 
lich eine Spitze S! Dies wird 
klar, wenn man sich p, durch 
Abwicklung der Evolute e von m 
erzeugt denkt, die bekanntlich 
eine zu m ähnliche Zykloide ist. 
Bei der Aufwicklung des ge- 
spannten Fadens f auf e beschreibt sein Endpunkt F das 
Profil p,, setzt einmal in-S mit einer Spitze senkrecht 
auf e auf und setzt dann die Bahn gegenläufig fort, bis 
er in P den Kreis p, oskulierend durchsetzt. Die Spitze S 
liegt mithin im Innern des Kreises p,, so daß dieser das 
Gegenprofil p, außer in P noch in einem Punkt @ trifft, 
wo p, dann eine stumpfe Ecke erhält (Bild 2). Die Zahn- 
Bild 1. lücke wird daher nicht von einem Halbkreis, sondern von 


1) 2. B. L. Burmester: „Lehrbuch der Kinematik“ Bd.TI, Leipzig 1888, S. 191f.; A. Schiebel: „Zahnräder“ 
I. T., Berlin 1922, S. 103. 
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einem etwas kürzeren Bogen Q9 gebildet, und vom eigentlichen Zahnprofil geht das Stück QSP 
verloren. Die Triebstockverzahnung läßt sich somit überhaupt nicht einwandfrei ausführen, 
wenn die Triebstockmitten auf dem Teilkreis liegen. 

Dieser Umstand ist es wohl, der Hartmann’) veranlaßt, die Triebstockverzahnung 
kurzweg und vollständig abzulehnen, „weil sie in der Nähe des Teilkreises unrichtig wird“. 

Nun, so schlimm ist die Sachlage wiederum nicht. Die Differenzen sind überaus klein, 
für Präzisionsübertragungen wird die Triebstockverzahnung sowieso nicht verwendet und ein 
gewisses Flankenspiel ist ja immer vorhanden. Überdies können — ja sollen sogar — mehrere 
Triebstöcke gleichzeitig im Eingriff stehen. 


8 2. Abschätzung des Flankenspiels. 


Das größte durch die Unstimmigkeit der Triebstockverzahnung verursachte fehlerhafte 
Flankenspiel ist dem Betrage nach gleich dem Abstand der Profilecke @ vom Bogen PS 
(Bild 2). Wir gelangen zu einer Abschätzung, indem wir die Entfernung der Profilspitze S 
vom Triebstockumfang p, berechnen. Da dessen Halbmesser s als Bogen MS auf die Evolute e 
aufgewickelt erscheint, handelt es sich um die Bestimmung des Überschusses des Zykloiden- 


bogens MS= s über die zugehörige Sehne MS=o. 

Wir setzen zu diesem Zwecke die Zykloide e in Parameterdarstellung an: 
xz=asinbr+bsin 
y=acosbr+bcosar 


Bezeichnen r, und r, die Rollkreisradien der Zykloide m, so geht diese in ihre Evolute e durch 
ähnliche Verkleinerung im Verhältnis r,:(2r,+ r,) über und es ist 


(r,+r,)r, 


Die Bogenlänge s beträgt 


das Quadrat der Sehne o 


®=2’+(y—a+b) | 
—=2a(a +b)(1—cosdbr) 


Da der gesuchte Unterschied s — o gegenüber s sehr klein ist, können wir mit ausreichender 
‚Näherung 


(4). 


s’ — 0? 


setzen und erhalten so 


8 (a — b)? 


Entwickeln wir schließlich (6) in eine Potenzreihe nach ı und auf Grund von (3) r in eine 
Potenzreihe von s, so erhalten wir das gesuchte Spiel als Reihe in s: 


a+b)? 2a? — 5ab+2b 


Das erste Glied dieser Entwicklung läßt sich einfacher auch so ableiten, daß man die 
Zykloide e durch ihren Scheitelkrümmungskreis ersetzt (Radius 0,=4ab/(a+ b)) und den 
Unterschied zwischen dem Kreisbogen s und der zugehörigen Sehne berechnet. — 

Drücken wir endlich a und b nach (2) durch die Wälzkreisradien r,, r, aus, so haben 
wir in erster Näherung: 


s—0= 


b(l-cosar) (6). 


(2r, +r,)'-s’ 


Für den in Bild 1 zugrundeliegenden Fall r,=2, r,=3, s=0,4 liefert das ein Flankenspiel 
von 8— o=0,00045, das sind rund 0,11 °/ des Triebstockhalbmessers. 


2) W. Hartmann: „Die Maschinengetriebe‘ Bd. I, Stuttgart 1913, S. 390. 
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Bild 3 zeigt eine nomographische Darstellung des Spiels auf Grund von Formel (8) in 
Gestalt einer Leitertafel von projektivem Typ. Sie ist in erster Linie für Übersetzungen 
n=r,[r, Z1 gedacht, da der Triebstockhalbmesser auf r, bezogen ist; für n<.1 empfiehlt 
sich der Bezug auf r,. Man gewinnt jedenfalls einen guten Überblick über die Größen- 
ordnung des fehlerhaften Spiels. 

Besondere Erwähnung verdienen noch die beiden Grenzfälle, in denen eines der Räder 
in eine Zahnstange übergegangen ist. Für r,—=x wird 

und die Leitertafel kann auch dafür verwendet werden (n=»). Für r,—», das ist der in 
der Praxis häufige Fall, in dem die Zahnstange die Triebstöcke trägt, gilt 
s’ 


Diese Abhängigkeit ist durch die Doppelleiter in Bild 3 unten dargestellt. 


[7 
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Bild 3. 


$ 3. Behebung des Fehlers. 


Wir werfen nun die Frage auf: Um welchen Betrag h müssen die Triebstöcke ins 
Innere des Teilkreises gerückt werden, um eine einwandfreie Verzahnung zu gewährleisten? 

Daß auf diese Weise eine Abstellung des in $ 1 beschriebenen Mißstandes möglich ist, 
kann man sich leicht klar machen: Verschiebt man die Triebstöcke um ein beträchtliches 
Stück h gegen die Rollkreismitte O,, so beschreiben ihre Mittelpunkte M bei der Rollung 
gestreckte Trochoiden, deren innere Parallelkurven im Abstand s keine Singularitäten 
aufweisen, wenigstens nicht an den kritischen Stellen. Hat man hingegen h zu klein ge- 
nommen, so löst sich die Oskulation bei P in zwei weitere Rückkehrpunkte auf und das 
Profil p, weist in dieser Umgebung nun sogar 3 Spitzen auf. Der Grenzfall — und gerade 
dieser muß uns interessieren — besteht nun darin, daß zwei dieser Spitzen sich vereinigen, 
und zwar zu einem sogenannten Gipfel- oder Spitzpunkt. . 

Wir können uns die Verhältnisse noch in der Weise anschaulich machen, daß wir die 
verschiedenen Parallelkurven einer gestreckten Trochoide m betrachten (Bild 4). Deren 
Evolute e besitzt eine Asymptote w, dem Wendepunkt von m entsprechend, und einen 
Rückkehrpunkt R, hervorgerufen von einem „Zwischenscheitel* der Trochoide. Ist mithin 
der Abstand s klein, so begleitet die Parallelkurve I die Trochoide auf der Innenseite ohne 
Singularitäten; ist s jedoch groß, so erhält die Parallelkurve III 2 Spitzen (man verfolge die 
Abwicklung der Evolute); im Grenzfall II geht sie ohne äußere Auffälligkeiten durch R, 
besitzt hier aber als höhere Singularität einen Gipfelpunkt. 

Diese Überlegung bietet auch schon die Handhabe zur rechnerischen Beantwortung der 
eingangs aufgeworfenen Frage. Wir setzen die Trochoide in Parameterdarstellung an: 

z=asinbr— csinar 


(9). 


y=acosbr—ccosar 
Hierin it a=r,+r,, b=r, und e=r,—h der Zentralabstand 0, M des Triebstockes. 
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Wir suchen nun jene Punkte, deren Krümmungshalbmesser o dem Betrage nach gleich 
dem Triebstockradius s ist; diese führen nämlich zu den Spitzen der Parallelkurve (vgl. III 
in Bild 4). Es ist 


_ —2be-cos(a — br]? 
(ac +b) (10) 


und die Bedingung o= + s führt bei der Bestimmung des Parameters r auf die kubische‘ 


Gleichung 

in der abkürzend 


gesetzt ist. 


Wir sehen, die Ermittlung der Spitzen für die Parallelkurve einer Trochoide ist im 
wesentlichen eine kubische Aufgabe. Sollen zwei Spitzen zusammenrücken, so muß (11) 
eine Doppelwurzel haben, d.h. es muß neben (11) auch noch die durch Ableitung nach t 
entstehende Gleichung 


bestehen. Entfernung von t liefert uns dann die Diskriminantengleichung | 
Über (12) ergibt sich daraus unsere gesuchte Grenzbedingung: 


(15). 


Damit ist die Formel für den Mindestbetrag h gefunden, um den die Triebstöcke gegen 
die Radmitte versetzt werden müssen, so daß das Auftreten von Spitzen beim Gegenprofil p, 
mit Sicherheit verhütet wird. Das Gegenprofil enthält nun keinerlei Kreisbogen mehr, 
sondern besteht zur Gänze aus einer Parallelkurve 
2% der Mitteltrochoide m. 
Da h sehr klein ist, kann das Glied Ah? in 
Formel (15) vernachlässigt werden. 
17% Für die in Bild 1 zugrundegelegte Annalıme 
r,=2, r,=3 s=0,4 ergibt sich h = 0,0068, d. i. 
0,34 °/s vom Teilkreishalbmesser r.. 
[93 Bild 5 gibt eine nomographische Darstellung 
der Funktion (15) in Form einer Leitertafel; 
Triebstockhalbmesser s und Verschiebung h sind 
dabei auf den Wälzkreisradius r, bezogen. — 
1005 Erwähnt werden müssen wieder die beiden Grenz- 
fälle ,=® und r,=x. Ersterer wird noch durch 
die Leitertafel erfaßt (n=x), letzterer führt auf 


\ 


gr 


02 + 


Verschiebung h in % von T, 


03+.0” 


Verschiebung in %ovon, 48? 
7 02 03 2 
und wird durch die Doppelleiter in Bild 5 unten 
Bild 5. dargestellt. 


Abschließend kann also festgestellt werden: Der Fehler der Triebstockverzahnung ist 
außerordentlich klein. Das durch ihn bedingte Flankenspiel liegt in der Größenordnung 
des praktisch stets vorhandenen Spiels. Durch geringfügige Verschiebung der Triebstöcke 
gegen die Radmitte läßt sich der Fehler vollständig ee ohne dafi die Eingriffs- 
verhältnisse merklich verschlechtert würden. 479 
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Lösungen der inhomogenen Mathieuschen Differential- 
gleichung mit periodischer Störfunktion beliebiger Frequenz. 
(mit besonderer Berücksichtigung der Resonanzlösungen). 

Von Gertrud Kotowski in Berlin. 
(Mitteilung aus dem Vierjahresplan-Institut für Schwingungsforschung, Berlin.) 


Die vorliegende Arbeit ist der systematischen und vollständigen Unter- 
suchung der bis jetzt nur in allgemeinen Zügen bekannten Lösung der 
inhomogenen Mathieuschen Differentialgleichung mit und ohne Dämpfungs- 
glied gewidmet. Es wird vor allem nach der Form der Lösungen (mit 
Einschluß der Resonanzlösungen) und nach der Lage der Kurven der 
Resonanzstellen gefragt. Da sich unendlich viele Kurven der Resonanz- 
stellen ergeben, wird, bei gegebener Störfunktion und Dämpfungsfreiheit, 
eine numerische Untersuchung für gewisse Werte der Parameter vorge- 
nommen, um aus der „Breite“ der Resonanzbereiche eine Abschätzung für 
die Bedeutung der Resonanzstellen mit wachsender Ordnungszahl zu ge- 
winnen. Ferner werden für einige Werte des Dämpfungskoeffizienten 
„Dämpfungskurven“, die zugleich die Kurven der Resonanzstellen bei ge- 
wissen Störfrequenzen sind, berechnet. 


A. Die inhomogene Mathieusche Differentialgleichung ohne Dämpfungsglied. 


I. Lösung im Inneren der Gebiete (u+"5*) der Struttschen Karte. 


a) Herleitung der Form der erzwungenen (partikularen) Lösung. Die 
inhomogene Mathieusche Differentialgleichung mit periodischer Störfunktion lautet 


wobei f(p,t) eine periodische Funktion von £ mit der Kreisfrequenz p sein soll. Setzt man 
wie bei der homogenen Differentialgleichung 


p 
pi=g8, 
so erhält die Gl. (1) die Form (wenn Striche Ableitungen nach x bedeuten) 
mit 
a? 


Die Mathieusche Differentialgleichung ist eine gewöhnliche, lineare Differential- 
gleichung zweiter Ordnung. Es gilt der Satz, daß, wie auch immer die Koeffizientenfunktionen 
einer solchen Differentialgleichung aussehen, man das allgemeine Integral der inhomogenen 
Gleichung bei Kenntnis eines Fundamentalsystems der homogenen Differentialgleichung stets 
bestimmen kann. Die allgemeinste Methode zu dieser Bestimmung ist unter dem Namen der 
„Variation der Konstanten“ bekannt. Sie erlaubt die Bestimmung der gesuchten Lösung mit 
Hilfe einer einfachen Quadratur. 

Seien y, (x) und y, (x) ein Fundamentalsystem der zu (2) gehörigen homogenen Gleichung, 
dann bestimmt sich das allgemeine Integral y, (x) der inhomogenen Gl. (2) aus der Formel 


x x 


mit 

I=Y — 9, Yı const =F0. 
Hierbei sind die unteren Grenzen der Integrale beliebig; sie entsprechen den willkürlichen 
Integrationskonstanten in der Lösung der homogenen Differentialgleichung. Es interessiert 


aber nur der Anteil der Lösung, der von der oberen (variablen) Grenze herrührt. Er soll im 
folgenden stets als Y,arı bezeichnet werden. 
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Da wir es mit einer linearen Differentialgleichung zu tun haben, dürfen wir als Stör- 
funktion in der zu behandelnden Differentialgleichung (2) eine einzige Harmonische 


betrachten. 
Die Fundamentalintegrale der zugehörigen homogenen Differentialgleichung besitzen (so- 


lange man sich nicht auf einer Grenzkurve der Struttschen Karte [l1, S. 216]'), Bild 1, be- 
findet) die Formen 


La, e'’* 
und 


+» 


P | / Z Bild ı Stabile 
d labile Lö 
Di ialglei- 
N 7) | 
2) 
Z / 


/ % 
G 
/ 
Bild 2 (rechts). Ver- 3 3 
halten des charakte- 
ristischen Exponenten 
“=u(l,y) auf einer 7 
Parallelenzur 4-Achse. 


(Hierin ist «= u(4,y), der charakteristische Exponent [1, S. 216], für die Zwecke der vor- 
liegenden Untersuchung eindeutig so festgelegt, wie es Bild 2 zeigt.) Führt man die Aus- 
drücke (4) für die Störfunktion und (5) für die Fyndamentalintegrale in die Formel (3) ein 
und vertauscht hierin Summen- und Integralzeichen, was bei der Art der Reihenentwicklungen 
der Lösungen der homogenen Differentialgleicnung erlaubt ist, so ergibt sich die Lösung 
Ypart in der Form 


x 
iv+i— 
(6). 
und — u): 
Nach Ausführung der Integration erhalten wir daraus, wenn für alle » der Ausdruck 
ist, 
+® +® 
Pe (8). 
-» 


1) Die eckigen Klammern im Text weisen auf das Schrifttumverzeichnis am Schluß der Arbeit hin. 
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Multiplizieren wir aus, so ergibt sich %part in der Form 


ivx 


da e*.e”“*—=1 und das Produkt zweier periodischer Funktionen derselben Frequenz im 
allgemeinen Falle wieder eine periodische Funktion dieser Frequenz ist. 


Das Ergebnis läßt sich also in der folgenden Weise aussprechen: Ist. (i» +i 2, + „)+0, 


so erhalten wir im allgemeinen Falle als Form des partikularen Integrals der inhomogenen 
Mathieuschen Differentialgleichung im Innern der Gebiete der Struttschen Karte (wie 
bei der Lösung der homogenen Gleichung) eine periodische Funktion zweiter Art; hier ist je- 
doch auch der erste Faktor stets periodisch und daher stets beschränkt; er hat die Frequenz 
der Störfunktion, während der zweite Faktor wieder die Frequenz der Belegungsfunktion 
aufweist. 

“Die Form des partikularen Integrals (9) legt es nahe, zwei wesentliche Fälle der Größe 


L (Störfrequenz) zu unterscheiden: 


und (n=ganze Zahl einschließlich der Null). 
Der Fall 2 wird im folgenden stets nochmals unterteilt in 


gr 5=75 


Wir unterscheiden also bei der weiteren Untersuchung stets die drei Arten von Stör- 


frequenzen: Die Störfunktion hat 


a) irgendeine, von n und +! verschiedene „Kreisfrequenz“, 


ß) die „Kreisfrequenz“* n, 
2 


y) die „Kreisfrequenz* 
(n= ganze Zahl einschließlich der Null). 


Wir wollen zunächst die Form der Lösungen im Innern der Gebiete in den drei 
Fällen a), 8) und y) gesondert angeben unter der Voraussetzung, daß die Bedingung (7) er- 
füllt ist. 

Im Fall a), also wenn die Kreisfrequenz der Störfunktion irgendeine beliebige reelle 


von n und et verschiedene Zahl ist, ist das partikulare Integral eine beschränkte perio- 
dische Funktion zweiter Art, deren erster Faktor die Frequenz der Störfunktion und deren 
zweiter Faktor die Frequenz der Belegungsfunktion hat (ist das Verhältnis von 5 zu n 


rational, so ist die Gesamtfunktion bekanntlich wieder periodisch). In dem zweiten periodi- 
schen Faktor werden im allgemeinen sämtliche Harmonische auftreten, auch wenn die Stör- 
funktion nur eine einzige Harmonische enthält. 

Im Fall 3), wenn also die Frequenz der Störfunktion gleich der der Belegungsfunktion 


oder gleich einem ganzzahligen Vielfachen davon ist, ergibt sich mit 5 =n für das partikulare 


Integral (9) die Form z 


Die partikulare Lösung ist also eine beschränkte, periodische Funktion mit der Frequenz der 
Belegungsfunktion. Selbst wenn die Störfunktion nur eine einzige Harmonische enthält (es 
ist auch »=0, also der Fall einer konstanten Störfunktion, eingeschlossen), enthält die er- 
zwungene Lösung im allgemeinen Fall sämtliche Harmonische mit Einschluß der Konstanten. 

Entsprechende Resultate ergeben sich im Fall y). Setzt man in der Form (9) des 


partikularen Integrals die „Kreisfrequenz“ Zn, ist also die Frequenz. der Störfunktion 


gleich der halben Frequenz der Belegungsfunktion oder gleich einem ungeraden ganzzahligen 


tör- 
(so- 
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Vielfachen der halben Frequenz, so erhält man als Lösung eine beschränkte, periodische 
Funktion mit der halben Frequenz der Belegungsfunktion. (Ganzzahlige Harmonische können 
in der Lösung nicht auftreten, sie enthält also sicher keine Konstante.) 

b) Die partikulare Lösung in den Fällen £) und y) ist gerade oder unge- 
rade, je nachdem die Störfunktion gerade oder ungerade ist. Wir wollen in 
diesem Abschnitt das bereits in der Überschrift angegebene Resultat kurz ableiten, und zwar 
nur für die Fälle 3) und y), für die die Lösung %,art eine rein periodische Funktion ist. Im 
Falle «), wo %yart im allgemeinen die Form einer beschränkten periodischen Funktion zweiter 
Art hat, ist die Unterscheidung in gerade und ungerade Funktionen zweiter Art (auch im 
Hinblick auf die später erfolgenden Untersuchungen) ohne wesentliches Interesse. 


Im Fall 5) (4=) erhielten wir die partikulare Lösung in der Form (10). Schreibt 


man statt der e-Funktionen die gleichwertigen Kreisfunktionen an, so ergibt sie sich in der 
Gestalt 


1 
Führt man diesen Ansatz in die inhomogene Mathieusche Differentialgleichung 


y= agcosox+ sing x 
0 1 


ein, so erhält man nach Durchführung des Koeffizientenvergleiches ein unendliches, inhomogenes 
System algebraischer Gleichungen. das in zwei zerfällt, eines zur Bestimmung der A, und 
ein zweites zur Bestimmung der B,. Die linken Seiten dieser Systeme stimmen mit denen 
zur Bestimmung der Fourierkoeffizienten der Mathieuschen Funktionen erster Art, mit 
dem Index » überein [2, S. 30]. Also entsprechen auch die Systemdeterminanten Jc„ und Ag, 
dieser homogenen Systeme den Nennerdeterminanten der inhomogenen Systeme zur Be- 
stimmung der A, und B,. Ist nun die Störfunktion gerade, sind also alle &, gleich Null, so 
wird das System zur Bestimmung der Koeffizienten B, homogen, die Determinante muß also 
zu Null werden, wenn eine nicht-triviale Lösung existieren soll. Da diese Determinante aber 
gleich der Determinante des Systems zur Bestimmung der Fourierkoeffizienten der Ma- 
thieuschen Funktionen erster Art, Ag„, ist, von der wir wissen, daß sie nur für die Para- 
meter der S,-Kurven zu Null wird [2, S. 41], so werden also sämtliche B, zu Null; bei ge- 
rader Störfunktion ist also die partikulare Lösung im Innern der Gebiete gleichfalls eine 
gerade Funktion. Der gleiche Gedankengang läßt sich für eine ungerade Störfunktion durch- 
führen. Die partikulare Lösung ist also bei ungerader Störfunktion überall im Innern der 
Gebiete gleichfalls eine ungerade Funktion. 

Entsprechende Resultate gelten im Fall y). 

c) Resonanzbedingung und Kurven der Resonanzstellen. Nur im Fall a) 
liegen die Kurven der Resonanzstellen im Innern der (stabilen) Gebiete. 

Die in den Abschnitten a und b abgeleiteten Formen der Integrale bleiben gültig, so- 
lange die Bedingung (7) nicht verletzt ist. Es ist nun als nächstes die Frage zu beantworten, 
ob und für welche Parameterwerte A und y die Bedingung (7) nicht mehr gilt, wann also 


ist. Wir wollen die Bedingung (12) im folgenden stets als „Resonanzbedingung“ bezeichnen. 
Für y=0 geht sie in die bekannte Bedingung für die lineare Differentialgleichung zweiter 
Ordnung mit konstanten Koeffizienten ohne Dämpfungsglied, über, da 


ist. Wie Formel (8) zeigt, wird die partikulare Lösung für jene Werte von u(A,y) (und »), 
die der Bedingung (12) genügen, unbeschränkt, da dann einige der Nenner zu Null werden. 
Alle jene zusammengehörigen Parameterwerte A und y, die zu einem Wert u (A,y) = constans 
gehören, für den die Bedingung (12) erfüllt ist, wollen wir als „Kurven der Resonanzstellen“ 
bezeichnen. (Bei einem festen Wert von v erhalten wir eine Kurve der Resonanzstellen 
für einen bestimmten Wert p/2=constans; da v eine ganze Zahl ist, die sämtliche Werte 
zwischen —o und + annehmen kann, erhalten wir für einen festen Wert p/2 = constans, 
also unendlich viele Kurven der Resonanzstellen.) Für die Wertepaare A und y der Kurven 
der Resonanzstellen wird die partikulare Lösung also unbeschränkt. Aus der Bedingung (12) 
bestimmen wir nunmehr die Lage der Kurven der Resonanzstellen in den Fällen «), £) und y). 

Wir wenden uns zunächst dem Fall.«a) zu. Da » eine ganze Zahl ist, & aber weder 


ganz- noch halbzahlig sein sollte, so kann Gl. (12) nur dann erfüllt sein, wenn auch . weder 
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ganz- noch halbzalılig, aber rein imaginär ist. Das kann aber nur für solche Wertepaare A 
und y der Fall sein, die innerhalb der stabilen Gebiete der Struttschen Karte 
[1, S. 216] liegen und für die 


ist. Auf den Kurven solcher zusammengehöriger Werte Ä und y innerhalb der stabilen Ge- 
biete der Struttschen Karte, die durch längs der Kurven konstante, der Gl. (13) genügende 
Parameterwerte « charakterisiert werden, wird die 


N partikulare Lösung also unbeschränkt; Gl. (13) stellt 
PA N j damit die Definitionsgleichung der Kurven der Re- 
! N (7 NN sonanzstellen im Falle «) dar. Wie man aus Gl. (13) 
% N N N weiter ablesen kann (und wie wir oben schon er- 
oc [5% N ZEIIIN wähnt haben), gibt‘ es selbst bei einem fest ge- 
N RN gebenen Verhältnis von 5 unendlich viele Kurven 
N N \ ie jeden Wert zwischen — und ‘+ annimmt. 
INN N SI N In Bild 3 sind schematisch die ersten Kurven der 
| N IR N Resonanzstellen für das Beispiel ot ein 
N N N gezeichnet. 

N Im Fall nimmt die Bedingung (12) mit 

N N die Form 
N N ... . (4 
[| N B an. Da » und n ganze Zahlen sind, so ist Gl. (14) 
nur dann erfüllt, wenn auch « eine imaginäre ganze 
de Zahl ist. Dies ist aber nur auf den ganzzahligen 
Fall  Grenzkurven der Struttschen Karte der Fall. Sie 
(5#5) für Beispiel (5; =+35) stellen also die Kurven der Resonanzstellen des 


(schematisch). 


Falles $) dar, und zwar sind, auch bei festem L=n, 
sämtliche ganzzahligen Grenzkurven Kurven der Resonanzstellen, da » ja alle ganzen Zahlen 
von — © bis + © durchläuft. 


Im Fall y) erhält man mit 2-"t ! aus Gl. (12) nach denselben Überlegungen 


das Ergebnis, daß für jedes n alle halbzahligen Grenzkurven Kurven der Resonanzstellen 
sind. In den Fällen $) und y) ist die Lösung im Innern der Gebiete also stets beschränkt. 

d) Herleitung der Resonanzlösung für den Fall a). Da die Kurven der Re- 
sonanzstellen in die stabilen Gebiete fallen, also « sicher verschieden von n ist, gilt das 
Fundamentalsystem 


(1) 


auch für die Wertepaare A und y, die der Resonanzbedingung (13) genügen. 

Damit bleiben auch alle Formeln, die wir unter (a) zur Herleitung des partikularen 
Integrals benutzten, gültig, solange wir noch nicht integriert haben. Formel (8) dagegen, 
also nach Ausführung der Integrale, bleibt nur für jene Werte von u (und ») richtig, für die 
die Bedingung (7) erfüllt ist. Diese Glieder geben nach dem Ausintegrieren wieder eine be- 
schränkte, periodische Funktion zweiter Art, also ein Integral. der Form (9), während die 


Glieder, für die ‚=-$5 + z ist, nach dem Ausintegrieren liefern: 


+» 
—ıe Du » a,e bzw. xe m.» 
i 


15 
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Das gesamte partikulare Integral ist also von der Form: 


Die Resonanzlösungen im Falle a) gehen also mit x’ gegen Unendlich. 


I. Lösungen auf den Grenzkurven R =") der Struttschen Karte. 


a) Herleitung der Form der erzwungenen Lösung. Wenn « die Werte "5: 


annimmt, so ist das zugehörige Fundamentalsystem der homogenen Differentialgleichung nicht 
mehr von der Form (15). Wenn a=n:-# ist, hat es bekanntlich die folgende Gestalt [1, S. 216] 


La, e'* | 
9 ()=ux Su, | 
wenn — i ist, lautet es: 
2 
x 
()+h,(@)=ßxe ? b,e'’* (18). 


Wir schreiben in (18) uns das Integral auch in der Zwischenform an, wie es vor dem Herein- 


multiplizieren mit et“*=e aussehen würde, weil diese uns in den folgenden Her- 
leitungen nützlich sein wird. Das partikulare Integral %)art berechnen wir wieder mit Hilfe 
der Gl. (3) (Variation der Konstanten). 


Wir führen für y,(x) und 4, (x) die Fundamentalintegrale (17) (im Falle u= ni) bzw. 


die Fundamentalintegrale (18) (im Falle PL. : ) in die Gl. (3) ein, vertauschen Summen- 


und Integralzeichen miteinander und integrieren aus; dann ergibt sich im Falle a), da hier 
sicher keiner der Exponenten der e-Funktionen unter dem Integralzeichen verschwinden kann 


(",n =ganze Zahlen, 5+3), für u=n-i: 
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bzw. für 


+® 
a ivx ivx 
Ynart = ßxe > a,e "—+e > b,e 
2n+1 
> > a, pie 2 
p\ 


Wie man leicht sehen kann, heben sich die mit x multiplizierten Glieder in den Formeln 
(19) und (20) heraus. Damit erhalten wir für die gesuchte Lösung bei den Werten u=in 
eine Funktion der Art 
x 
(21). 


Das gleiche Resultat ergibt sich für die Werte u= Be i, wenn wir noch beachten, daß 
2n+l, 2n+l,, 
das Produkt e ? .e ? = 
Im Fall a), also für +7: ist daher auch für jene Wertepaare A und y, die zu den 


Grenzkurven der Struttschen Karte gehören, die partikulare Lösung %»art beschränkt und 
von der Form (9), ein Resultat, das sich damit deckt, daß für die Grenzkurven die Resonanz- 
bedingung (13) nicht erfüllt ist. 


Im Fall ?), also für 5=n, erhalten wir in gleicher Weise das Ergebnis, daß die 


Form (10) des partikularen Integrals auf den halbzahligen Grenzkurven erhalten 
bleibt. Anders dagegen verhält sich hier die Lösung auf den ganzzahligen Grenzkurven, für 
die ja die Resonanzbedingung (14) erfüllt ist. Die in Formel (19) erhaltenen Integrale werden 


für den Wert v„=— n ungültig, während sie für alle anderen Werte von » weiterhin gelten. 
Wir schreiben uns den Anteil des partikularen Integrals für »=—n heraus und integrieren 
ihn gesondert. Das ergibt 
a +» +» x 
> 
Ausintegriert, ausmultipliziert und entsprechend zusammengefaßt, erhalten wir daraus 
Die Gesamtlösung ist also von der Form 
ie», 


sie stellt die „Kesonanzlösung“ des Falles £) dar. 
Entsprechend erhalten wir im Fall y) das Resultat, daß auf den ganzzahligen Grenz- 
kurven der Struttschen Karte die Lösung %,arı von derselben Form wie in den Gebieten 


ist, während die Lösungen auf den halbzahligen Grenzkurven (ur ) vom Typ 


15° 


2 
216] 
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sind; sie sind die „Resonanzlösungen“ des Falles y). 

Im Fall 5) und y) gehen also die „Resonanzlösungen“ mit x? gegen Unendlich, während 
sie im Fall a) nur mit x' gegen Unendlich gehen. 

Wir fassen die bisher erhaltenen Resultate kurz in dem folgenden Satz zusammen: Das 
partikulare Integral der inhomogenen Mathieuschen Differentialgleichung ohne Dämpfungs- 
glied mit periodischer Störfunktion wird nur in den Fällen, in denen die Resonanzbedingung 


(iv+ + u(d, „)=0 erfüllt ist, unbeschränkt. Für alle anderen Wertepaare A und y, ganz 


gleich, ob sie in die Gebiete oder auf die Grenzkurven fallen, bleibt die Lösung beschränkt. 
Diese Aussage gilt für alle möglichen Werte der Frequenz der Störfunktion (Fall a), $) und y)). 


b) Trennung der Kurven der Resonanzstellen in den Fällen ) und y), 
je nachdem, ob die Störfunktion gerade oder ungerade ist. Wir wollen hier 
noch kurz auf eine Einschränkung zu den Ergebnissen des Abschnitts Ic in den Fällen £) 
und y) aufmerksam machen, die dann auftritt, wenn die Störfunktion entweder nur gerade 
oder nur ungerade ist. Wir betrachten als Beispiel die partikulare Lösung auf den C„-Kurven 
im Fall 8). Die Störfunktion habe eine der beiden Formen 


1. und 2 . . . 2 2 2.0. (MB). 
Das Fundamentalsystem, das auf den ©„-Kurven gilt, kann statt nach (17) auch in der Form 


Yı a,cosvx, 
em 
0 1 


angeschrieben werden. Führen wir dieses Fundamentalsystem und die Störfunktion a, cosnx 
in die Gl. (8) ein, so ergibt sich nach Durchführung der Integration als Form der Lösungen 


Ypart= e,cosvx+x e,cos»2 . .... 


Setzen wir dagegen wieder das Fundamentalsystem (27), aber die Störfunktion f„sinnz in 
Gl. (3) ein, so erhalten wir zunächst als Ergebnis der Integration auch einen mit x gehenden 
Anteil. Zieht man aber von dem so erhaltenen Ausdruck für %part die mit einem bestimmten 
konstanten Faktor multiplizierte Fundamentallösung 4, (x) aus (27) ab, was stets erlaubt ist, 
so entfällt der mit x gehende Anteil. Wir erhalten also ein partikulares Integral, das be- 
schränkt ist und die Form 


hat. Entsprechendes erhält man auf den ungeraden ganzfrequenten Grenzkurven bei einer» 
ganzfrequenten geraden Störfunktion. 

Dieselben Resultate gelten im Falle y). 

Es zeigt sich also, daß in den Fällen 8) und y) für eine gerade Störfunktion nur die 
„geraden“ Grenzkurven, C,„ und C,„+1, Kurven der Resonanzstellen sind und für eine un- 


2 

gerade Störfunktion nur die „ungeraden“ Grenzkurven, S,„ und S;„+,, ein Resultat, das auch 
2 

mit dem Verhalten der Nennerdeterminanten des Abschnittes Ib in Einklang steht. 


c) Numerische Untersuchung der drei ersten Resonanzbereiche be- 
züglich ihrer „Resonanzbreite* für eine Störfunktion von der Form 
a,+a,cosz-+a,cos2x. Es handelt sich im folgenden um Untersuchungen, die zum Fall ß) 
gehören, und zwar bei einer geraden Störfunktion. 

Da die erzwungene Lösung %Ypart bei einer Störfunktion, die aus nur einer Harmonischen 
besteht, schon unendlich viele Resonanzstellen aufweist, erhebt sich die Frage nach der 
relativen Wichtigkeit der einzelnen Resonanzbereiche. Als Maß für diese Wichtigkeit wollen 
wir die „Breite“ der Resonanzbereiche einführen und dabei als „Resonanzbreite“ den Frequenz- 
bereich definieren, in dem die von den Störgliedern mit der Amplitude 1 hervorgerufenen 
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Zahlentafelii. y=1. 


2 4,5) 2 A, f 4A, 6) 
Go 
— 1,229 — 0,9993 — 0,839 — 1,0007 — 0,8 — 0,1166 
— 0,8  — 1,8494 — 0,6 | — 2,3408 — 0,6 — 0,2559 
— 0,6 — 3,6173 —0,5 — 4,5033 — 0,45 — 0,8896 
— 0,5 — 6,5033 — 0,45 — 7,8698 — 0,42 — 1,5697 
— 0:45 — 10,9665 — 0,4 — 26,9292 — 0,4 — 3,0788 
—0,4 -— 36,1616 — 0,379 oo — 0,379 oo 
— 0,379 © — 0,35 20,9470 — 0,35 2,4226 
— 0,35 27,0671 — 03 7,8400 — 0,31 1,0424 
— 0,3 9,7333 — 02 3,6767 — 0,3 0,9174 
—02 4,1918 — 01 2,5244 — 02 0,4406 
0 1,8741 {) 2,0000 0 0,2518 
0,225 1,0073 0,1 1,7137 05 0,2018 
0,5 0,5994 0,2 1,5462 0,8 0,2643 
0,8 0,2029 0,3 1.4502 1,0 0,4051 
0,9181 0 0,4 1,4045 1,1 0,5927 
1,0 — 0,2025 0,5 1,4006 1,18 0,9866 
1,1 — 0,6363 0,6 14378 | 12 1,1915 
1,145 — 1,0042 0,7 1,5233 1,25 2,5372 
1,2 — 1,9149 0,8 1,6753 1,293 © 
1,25 — 4,7160 1,0 2,4047 1,35 — 1,8831 
1,27 — 9,2040 1,1 3,3999 1,4 — 0,9920 
1,293 oo 1,2 6,5957 1,5 — 0,5042 
1,32 8,8080 1,25 13,7908 2,0 — 0,1451 
1,35 4,3914 1,293 © 
1,5 1.4958 1,35 — 9,8570 
1.655 1,0002 15 — 2,4874 
2,0 0,6425 1,75 — 0,9901 
2,0 — 0,5699 


erzwungenen Amplituden größer sind als + 1. Allerdings müssen wir dabei bedenken, daß 
von den Resonanzbereichen nur der Teil von praktischem Interesse ist, der in die stabilen 
Gebiete der homogenen Lösung fällt. Im folgenden wird als Beispiel die Differentialgleichung 
mit der Störfunktion a,+a,cosx+.a,cos2x und dem Wert y=1 untersucht. Die Größe 4° 
durchläuft dabei den Wertebereich —1 bis + 2,0. Das unendliche inhomogene System alge- 
braischer Gleichungen, das man zur Bestimmung der Fourierkoeffizienten A, erhält, wenn 


man mit dem Ansatz Ypart= 8 A,cosvx (AIb) in die inhomogene Mathieusche Differential- 
D 


gleichung mit der Störfunktion a,+a,cosz-+a,cos2x eingeht und einen Koeffizientenvergleich 
durchführt, konnte, wie sich bei Durchführung der Rechnung zeigte, wegen der guten Konvergenz 
bereits nach den ersten vier bzw. fünf Gleichungen abgebrechen werden. Es betrug nämlich 
bei den von uns gewählten Werten von A die Abweichung zwischen den Werten, die aus vier 
Gleichungen, und denen, die aus fünf Gleichungen berechnet wurden, nie mehr als '/a°J. 
Ferner zeigte sich, daß bereits der Fourierkoeffizient A, der partikularen Lösung im Ver- 
hältnis zum Fourierkoeffizienten A, vernachlässigt werden kann (A, beträgt höchstens noch 
10°, von A,). Es wurden deshalb nur die Koeffizienten A,, A, und A, bestimmt. Um die 
Tabulierung (und Aufzeichnung) möglichst allgemein zu halten, wurden in den Zahlentafeln 1 


bis 3 nicht die Fourierkoeffizienten selbst angegeben, sondern gewisse Quotienten a, die 


den Charakter von „Vergrößerungsfunktionen“ haben. Sie werden im folgenden auch stets 
als solche bezeichnet werden. Diese Bezeichnung wurde in Analogie zum Sprachgebrauch 
bei der linearen Differentialgleichung zweiter Ordnung mit konstanten Koeffizienten gewählt. 


Die Größe Zr bedeutet die Amplitude der »-ten Harmonischen in der partikularen Lösung, 


die durch die u-te Harmonische in der Störfunktion hervorgerufen wird, geteilt durch die 
Amplitude eben dieser «-ten Harmonischen. Da die Amplituden der partikularen Lösung 


Y»art nur linear von den Störamplituden abhängen, sind die Funktionen pr reine Funktionen 
von A und y. 
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Zahlentafel 2. y=1. 
A, A, 6) | 4,6) 
a, | a, | a, 
— 0.633 — 1,0588 — 0,734 — 1,0003 —05 — 0,2517 
— 0,6 — 11704 | — 0,7 — 1,0727 | — 0,45 0,4003 
—05 — 2,2517 | — 0,6 — 14055 | — 0,405 — 1,0076 
— 0,45 —3939 | —05 — 2,3517 | —04 — 1.2315 
—04 — 13.4646 | — 0,85 — 3,5414 — 0.39 — 2.2623 
— 0,379 — 0,4 1077171 —039 
—0,35 10,4735 | — 0,379 — 0,370 2,9585 
— 0,3 3,9200 | — 0,35 73315 | —0855 | 1,0388 
— 0,2 1,8384 | — 0.3 23520 | — 0,85 0,8479 
0 1.0000 — 0,23 0,9966 — 0,3 0,2752 
0,5 0,7008 | — 02 0,7353 | —02 | 0,0881 
0,8 0,8377 — 0,2544 0 
0,92 1,0037 0 0 0,5 —0,1009 
1,0 1.2024 0,1 — 0,1714 0,8 — 0,2114 
1,1 1,7000 0,2 — 0,3092 1,0 0,4051 
1,2 3,2979 0,3 — 0,4351 1,1 — 0,6519 
1,25 6,8954 0,4 — 0,5618 18 | — 1,1641 
1,293 © 0,5 — 0,7003 1,2 | — 1,4298 
1,32 — 10,6266 0,6 — 0,8627 | 15 | —3,1715 
1,35 — 4,9285 0,7 1,0663 1293 | 
1,5 — 1,2437 0,8 1,3403 | 15 2,5421 
1,55 — 0,9682 1,0 — 2,4047 10 | 1,3888 
1,75 — 0,4951 1,1 — 3,7399 1,5 0,7563 
23,0 — 0,2850 1,2 — 79148 2,0 0.2902 
1,25 — 17.2385 | 
| 1,293 | | 
| 1,35 13,3070 | 
| 1,5 3,7311 | | 
| 2,0 | 1.1399 | 


Sie ergeben sich nach der Cramerschen Regel als die Quotienten zweier Determinanten, 
die von derselben Ordnung wie die Gleichungssysteme, also von der vierten oder fünften 
Ordnung, sind. Zunächst sind nach den Inceschen Rechnungen [3] die Nullstellen der Nenner- 
determinanten, also die Unbeschränktheitsstellen der Lösung, die nach Abschnitt Ic auf den 
C„-Kurven der Struttschen Karte liegen müssen, bekannt. Für die „Breite“ der Resonanz- 
stellen der einzelnen Vergrößerungsfunktionen ist nun vor allen Dingen die Lage der Null- 
stellen der Zählerdeterminanten wichtig. Sie wurden deshalb zunächst bestimint (siehe auch 
Zahlentafeln 1 bis 3). Für einzelne von ihnen lassen sich die Kurven der Resonanzstellen in 
der A-y-Ebene allgemein (für alle Werte des Parameters y) angeben, und zwar liegen die 


Nullstellen von 4 stets auf den S„-Kurven der Struttschen Karte, die kleinste Null- 
o 
stelle von Zi, a und 2 stets auf der y-Achse, also bei A=0, während die beiden klein- 
1 2 
sten Nullstellen von 4 stets durch die Wurzeln der quadratischen Gleichung 


= 


gegeben sind, ganz unabhängig von dem Grad der zur Berechnung der Vergrößerungsfunktionen 
gewählten Annäherungsdeterminante. 


Die erhaltenen Zahlenwerte sind, wie schon erwähnt, in den Zahlentafeln 1 bis 3 
niedergelegt: Zur Bezeichnung ist noch zu sagen, daß der obere Index 5 angibt, daß der 
Zahlenwert aus fünfreihigen Annäherungsdeterminanten berechnet wurde. 


Wir wenden uns nunmehr einer kurzen Diskussion der Resultate zu. Für y=1 fallen 
in den betrachteten Bereich zwei Resonanzstellen. An der ersten Resonanzstelle (Zahlen- 
‘ tafeln 1 bis 3) bei 4,=— 0,379 ist der Resonanzbereich nicht nur für die Vergrößerungs- 
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Zahlentafel3. y=1. 
4, | 4; | | a, | 
| | 

—05 | | — 0,4008 | — 015 | — 0712 | 1,8 — 2,3383 
| —0405 — 1,0076 | - 04 — 0,3695 | 1,298 
— 0,41 — 1,044 | —040 — 1,2315 | — 0.384 — 0,7660 | 1,3 3,4833 
| —1534 | —09 | | | | 1381 1.2043 
—.0,39 | — 0,881 | — 1,4002 | 1,83 0,3544 
— 0,379 © | — 2176 | 185 | 0,1086 
— 0,37 4,0217 | — 0,355 1,0388 | — 0,379 | 1,366 
— 0,35 12113 | — 0,35 0,8479 | — 0,378 5,7568 | 1,5 | — 0,2521 
— 0,345 1.0348 | — 03 0,2752 | - 0,377 12338! 20 | —04332 
0,4587 | — 02 0,0881 — 0,376 1,0605 | 
— 02 0,2203 0 v | — 0,375 0,4021 | 

N) 0,1259 0,5 — 0,1009 | — 0,372 0,2169 | | 

0,5 0,1009 0,8 — 0,2114 | — 0,370 0,1104 | 

0,8 0,1322 1,0 — 0,4051 | — 0368 0,0344 | 

1,0 0,2026 1,1 — 0,6519 — 0,366 0 | | 

11 0,2964 1.18 — 1,1641 | -- 0,364 — 0,0330 | 

1,2 0,5958 1.20 — 1,4298 | — 0,360 — 0,0770 | | 

1,24 1,0304 1,22 — 1,8407 | — 0,355 — 0,1110 | 

1,25 1,2686 1,24 — 2,5613 | — 0,350 — 0,1332 | | 

1,27 2,3520 1,25 — 3,1715 | —03 — 0,2018 | | 

1,28 3,9828 1.27 — 5,9722 | — 02 — 0,2291 | | 

1,285 6,6481 1,293 © 0 — 0,3518 | 

1,293 © 1,31 8,4008 05 - 0,3027 | 

1,3 — 7,9166 1,32 5,2985 0,8 — 0,3489 | 

1,31 — 3,2064 | 1,85 2,5421 1,0 — 0,4051 | 

1,32 - 2,0070 | 1,97 1,9006 1,1 — 0,4623 

1,35 —0916 | 15 0,7563 1,2 — 0,6196 | 

1.37- — 0,6936 | 2,0 0,2902 15 —0815 | 

15 — 03521 | 16 | — 1,1357 | 

2,0 — 0,0726 127 | — 1,4780 


funktion 4 sehr breit, was zu erwarten war, sondern auch für die Vergrößerungsfunktion a, 
o o 


Dasselbe Bild erhalten wir 


4, 


während er für die Vergrößerungsfunktion ” sehr schmal ist. 
o 


für die Vergrößerungsfunktionen &. doch sind hier die Resonanzbreiten etwas geringer ge- 


worden. Dagegen sind alle Resonanzbreiten der von der Störamplitude a, herrührenden Ver- 
größerungsfunktionen sehr klein. 


Bei der zweiten Resonanzstelle, A, = 1,29, haben sich die Rollen von Zu und je ver- 


tauscht ; hier sind die Resonanzbreiten der Vergrößerungsfunktionen Zu noch bemerkenswert 


groß. Wie zu erwarten, ist die Resonanz für die Vergrößerungsfunktionen a hier stärker 


u 

ausgeprägt als an der ersten Resonanzstelle, wie auch alle durch die Störamplitude a, her- 
vorgerufenen Vergrößerungsfunktionen hier größere Resonanzbreiten aufweisen als an der 
ersten Resonanzstelle. 


Für größere Werte von A zeigt das Verhalten der Lösung %,art praktisch nur noch ge- 
ringe Unterschiede gegenüber dem Verhalten der Lösung der entsprechenden Differential- 
gleichung mit konstanten Koeffizienten (also für y=0). Wir stellten beim Berechnen der 
obengenannten Vergrößerungsfunktionen im Intervall A=2,5 bis A=8 fest, daß die Resonanz- 


breite der Vergrößerungsfunktion n in der Umgebung der dritten Resonanzstelle bei A, = 4,0353 
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bereits bei weitem die aller anderen Ver- 
— größerungsfunktionen überwiegt. Die Zahlen- 
sa. werte für diesen Bereich von A wurden des- 
7 halb nicht mehr mit aufgenommen. 
a In den Bildern 4 und 5 ist der Ver- 
Örenzen des stabilen Bereiches: A 
-93785_ und -93#77 lauf der Vergrößerungsfunktionen für 
die beiden ersten Resonanzstellen aufge- 
N zeichnet (Zahlentafel 1 bis 3). Die schraf- 
Av, fierten Gebiete sind hier stabil; die strich- 
punktierten Linien geben die Resonanz- 
stellen an. Wir wählten gerade dieses 
ı Beispiel, um zu zeigen, wie stark die Re- 
A 
sonanz von = noch ausgeprägt ist. Hieran 
1 1 

läßt sich nämlich, von einem etwas anderen 
Gesichtspunkt aus als es üblich ist [4], die 
Bis 4. Auswanderungserscheinung des Pendels mit 
horizontal erschüttertem Aufhängepunkt ab- 
2 lesen, ein physikali- 
Resonanzstelle 1. u scher Vorgang, der, 
A für nicht zu große 
Winkel, durch eine 
inhomogene Ma- 
thieu sche Differen- 
tialgleichung mit der 
Störfunktion a, cosx 
(mit der ursprüng- 
707 _Mittellage) beschrie- 

ben wird. (Für kleine 
Werte von A und y 
mit A<y, also in 
jenem Wertebereich 
von Aundy, für den 


bis jetzt untersucht 


Bild 5. wurde, ist der Unter- 


schied zwischen 4 und a im Bereich der ersten Resonanzstelle noch viel stärker aus- 
1 1 
geprägt, da 4 mit y, A mit A linear geht.) 
1 1 
B. Einfluß der Dämpfung. 
I. Lösung im Innern der Gebiete. Resonanzbedingung. Resonanzlösungen. 


Wir wollen im folgenden untersuchen, ob der Einfluß der Dämpfung auf die partikulare 
Lösung der inhomogenen Mathieuschen Differentialgleichung ein ähnlicher ist wie bei der 
Differentialgleichung mit konstanten Koeffizienten, d. h. vor alien Dingen, ob „scharfe“ Re- 
sonanzspitzen, also Unendlichkeitsstellen der Lösung, hier gleichfalls unmöglich sind. 


Wir schreiben die inhomogene Mathieusche Differentialgleichung zunächst in der Form 


diese Differentialgleichung kann bekanntlich durch die Transformation y=ze" °* in die 
folgende Differentialgleichung für z verwandelt werden: 


mit 
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Wir wollen die Lösung 2,.rt zunächst in den Gebieten bestimmen; das Fundamental- 
system ist also nach (5) gegeben. Berechnen wir die Lösung 2,.rt von Differentialgleichung (Ob) 
wie stets in der vorhergehenden Untersuchung über die Variation der Konstanten, so ergibt 
sich nach dem Ausintegrieren: 


— utiv 
a1), 
+® +» 
solange nur die Bedingung 
. 


für alle « und » erfüllt ist. Multiplizieren wir aus und gehen wieder zur ursprünglichen 
Koordinate über, so erhalten wir die Lösung in der Form 


+» 


ivx 
Yvart =de, 


Es ergibt sich also hier genau wie bei Abwesenheit von Dämpfung für die Lösung Yparı im 


Fall a) (#3) eine beschränkte periodische Funktion zweiter Art, in den Fällen ß) 


und y) eine beschränkte periodische Funktion mit der ganzen bzw. halben Frequenz der 
Belegungsfunktion, solange nur Bedingung (32) stets erfüllt ist. 

Wann ist nun diese Bedingung nicht mehr erfüllt? Die Ungleichung (32) erhält reelle 
und imaginäre Glieder. Damit sie nicht mehr gilt, müssen die beiden Gleichungen 


und 


einzeln erfüllt sein. Aus diesen beiden Resonanzbedingungen kann man sofort die folgenden 
Ergebnisse ablesen: In den stabilen Gebieten der Differentialgleichung (30b) kann Be- 
dingung (35) nie erfüllt sein, da hier « rein imaginär ist. Die Lösung %)art ist also in den 
stabilen Gebieten in allen drei Fällen «), $) und y) besehränkt. Die Kurven der Resonanz- 
stellen. in den stabilen Gebieten im Fall a), die wir für Dämpfungsfreiheit erhielten, ver- 
schwinden also; die Dämpfung übt hier einen ähnlichen Einfluß aus wie bei der ent- 
sprechenden Differentialgleichung mit konstanten Koeffizienten. In den instabilen Gebieten 
der Differentialgleichung (30b) ist im Fall a) die Bedingung (34) nicht mehr erfüllbar, die 
Lösung Ypart also überall in den instabilen Gebieten beschränkt. 

Im Falle £) lassen sich die Gl. (34) und (35) für jene Kurven in den „ganzzahligen“ 
(in Bild 1 mit 1 bezeichneten) Instabilitätsintervallen erfüllen, für die ö_ gleich dem reellen An- 
teil von «, den wir mit dem Buchstaben a bezeichnen wollen, ist; in den halbzahligen (in Bild 1 
mit '/, bezeichneten) Instabilitätsintervallen dagegen nieht, da der Imaginärteil von . hier halb- 
zahlig ist. Die Kurven a=Ö in den instabilen Gebieten der Lösung in z sind zugleich die 
Grenzkurven der homogenen Lösung der Differentialgleichung (30a), also der Lösung in y, 
wie man leicht verifizieren kann, wenn man sich die homogene Lösung der Differential- 
gleichung (30b) hinschreibt und dann mittels der Transformation y=2:e”°* zur Koordinate y 
übergeht. 

Für jene Werte von » und a, für die die Resonanzbedingungen erfüllt sind, ist aber 


Formel (31) nicht mehr gültig. (Man setze in der Formel I =n). Wir bestimmen diesen 


Anteil des Integrals gesondert und erhalten ihn in der Form (für »=—n+ Imu=m und 
a=6; Im u=ganze Zahl) 


+» 
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In der ursprünglichen Koordinate y geschrieben, ergibt sich also die partikulare Lösung Yyarı zu 


+® 


Entsprechend läßt sich im Falle y), also für I nd ‚ zeigen, daß die partikulare Lösung 
auf jenen Kurven in den halbzahligen instabilen Gebieten, für die a=ö ist, die Form hat 


2 2 
Yvart = +2 d, „+1 (38). 
2 2 


Die Lösungen (37) und (38) stellen die Resonanzlösungen des Falles $) und y) dar. , 

Wir erhalten also für die inhomogene Mathieusche Differentialgleichung mit Dämpfungs- 
glied das Ergebnis, daß, im Gegensatz zu den Lösungen der Differentialgleichung mit kon- 
stanten Koeffizienten, auch bei Anwesenheit von Dämpfung „scharfe“ Resonanz: 
stellen der Lösung (Unendlichkeitsstellen) auftreten. Und zwar ergeben sie sich in. den 


Fällen 8) und y), 5-7 also bei Störfrequenzen, die bei Abwesenheit von Dämpfung zu 
solchen Resonanzlösungen führen, die mit x? gegen Unendlich gehen. 


I. Lösungen auf den Grenzkurven a=0. 


Wir wollen jetzt noch kurz angeben, wie sich die Lösungen, die zu den Wertepaaren Ä 
und y der Grenzkurven (Grenzkurven für die Lösungen: in z) gehören, verhalten. Wir können 
hier ohne wesentliche Veränderungen die Überlegungen aus Abschnitt A IIa übernehmen, 
wenn wir uns stets statt der Koordinate 4 die Koordinate g und statt der Störfunktion 


iP x +9) x 

e die Störfunktion e' 2 gesetzt denken. Führen wir die Rechnungen hierfür. durch, 
so sehen wir, daß die Glieder mit x in den Ausdrücken für 2yart sich stets gegenseitig auf- 
heben, und zwar in allen.drei Fällen a), 8) und.;y) sowohl auf den ganz- als auch auf den 
halbzahligen Grenzkurven (a=0); dies deckt sic# Mit der Tatsache, daß in allen drei Fällen 
die Resonanzbedingungen (34) und (35) für die Parameter der Grenzkurven nie zusammen 
erfüllt sind. Die Form der Lösungen ist wieder durch Formel 33) gegeben; in den Fällen ß) 


und y) wird die Lösung daraus durch Einsetzen von 5=n bzw. nt s gewonnen. . 


II. Berechnung und Aufzeichnung‘ von „Dämpfungskurven“ in den drei ersten instabilen 
Gebieten nullter, halber und erster Ordnung für die Werte a?=ö?=10"*; 10°°?; 10°°; 
5-10”? und 10°'. 


Obwohl der Einfluß der Dämpfung auf die Veränderung der Stabilitätsgebiete und da- 
mit auf die Verschiebung der Grenzkurven zwischen stabilen und instabilen Gebieten qualitativ 
längst bekannt ist [l, S. 217], so findet man trotzdem in der vorliegenden Literatur nur ver- 
streut und unsystematisch Zalılenmaterial dazu [5]. Der Grund ist wohl in den großen Rechen- 
schwierigkeiten zu suchen, mit denen die Bestimmung der Kurven konstanter Werte von 
Reu(i,y)=a=constans in den instabilen Gebieten verbunden ist. Da diese Kurven 
Re u=a= constans stets bei einer ganz bestimmten Größe der Dämpfungskoeffizienten, näm- 
lich für a=ö, als neue Grenzkurven ig Erscheinung treten, so haben wir sie im folgenden 
auch als „Dämpfungskurven“ esciche Weil die Bedeytung dieser Kurven als Grenzkurven 
der Stabilitätsgebiete und zugleich als „Kurven der Resonanzstellen“ (siehe Teil BI) bei ge- 
gebener Dämpfung klar auf der Hand liegt, haben wir für einige Werte von a?—=ö?, nämlich 
für 10°; 5-10? und 10', diese Kurven bis zum Wert y=1 im in- 
stabilen Gebiet halber Ordnung und y=2,5 im instabilen Gebiet erster Ordnung bestimmt. 
Im instabilen Gebiet nullter Ordnung (siehe Bilder 1 und 2; der konstante Imaginärteil 
von « ist hier null) wurde wegen der verhältnismäßigen Kleinheit der von uns gewählten 
Werte a?=? nur die Kurve für a =ö’=10"' bis zum Wert y=1 berechnet. 

Die Kurven wurden im Gebiet halber Ordnung aus vierreihigen, im Gebiet nullter und 
erster Ordnung aus fünfreihigen Näherungsdeterminanten der Hillschen Determinante nach 
dem Hornerschen Näherungsverfahren bestimmt. Für das Gebiet halber Ordnung liefert 
bereits die quadratische Näherung mit y„in =2a=2Öd für die von uns gewählten Werte 
von a eine sehr gute Näherung. Es ergeben sich auf Grund dieser Rechnungen die in den 
Zahlentafeln 4 bis 6 angegebenen Werte und die in Bild 6 aufgezeichneten Kurven. Die Ge- 
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Zahlentafel 4. Instabiler Bereich 1/2. Ordnung. 
T 
a? = — 10-t a? = 5? — 10-3 | a — — 10-2 
0,0201 | 02488 | 0,250 | 0,0083 | 0,2477 0,2495 | 0,1998 0,2342 | 0,8358 
0,0203 0,2482 | 02518 | 0,0635 | 0,4676 0,2507 | 0,20 0,2332 | 0,2393 
0,0211 0,2469 | 0,2521 | 0,064 | 0,2437 0,570 | 021 | 0,2248 | 0,2452 
0,0218 0,245 | 0,2540 | 0,10 0,2120 0,2830 | 0202 | 0,2208 0,2495 
0,05 0,2270 | 0,2780 0,30 | 0,0920 0,3840 | 0,205 | 0,2124 0,2575 
0,10 0,2040 | 0,2970 | 0,50 | 0,0280 0,4640 0,30 ' 0,1180 | 0,3390 
0,30 0,0900 0,3870 | 10 | — 0,3460 0,5970 0,50 0,0150 | 0,4390 
050, —0,0850 | 0,4650 1,0 0,3450 | 0,5840 
100 1 0,3460 | 0,5980 | | 
y | y A, 
0,4614 0,1738 0.1760 0,6750 0,0994 0,1010 
0,4615 0,1708 0,1792 0,70 0,0123 0,1801 
0,4625 0,1595 0,1868 0,75 0,0537 0,2356 
0,4650 0,1548 0,2017 1,0 — 0,2940 0,3720 
0,470 0,1336 0,2118 | 
0,50 — 0,0790 0,2620 | | 
1,0 — 0,3270 0,4940 | | 
Zahlentafel 5. Instabiler Bereich erster Ordnung. 
—10-* | 10- | a2 = = 
| 
1 | | Y | 2 4, | 
0,2855 | 1.0119 | 1,0138 0,5168 1,0363 1,0379 | 0,96973 1,0881 1,08995 
0,2860 1,0110 1,0147 0,5170 1,0346 1,0395  0,96975 1,0873 1,0908 
0,2870 1,0099 | 1,0160 0,5180 1,0315 | 1.0429 0,9698 1,0860 1,0920 
0,2900 1,0082 | 1.0183 | 0,520 1.0281 1,0468 0,970 1,0834 1,0947 
0,30 ' 1,009 | 1,0233 | 0,580 1,0197 1,0579 0,990 1,0432 1,1405 
0,35 0,9973 1,0404 | 0,55 1,0104 | 1,0723 1,0 1,0335 1,1530 
0,40 | 0,9920 | 1,0562 | 0,60 0,9962 | 1,0999 | 1,1 0,9776 1,2365 
0,50 0,9825 | 1,0894 | 0,80 0,9598 | 1,199 | 12 | 0,9411 1,2996 
0,60 0,975 | 1,1260 0,85 0519 | 12152 | 15 | 0,8527 1,4590 
0,80 0,9485 1,2060 0,90 09433 | 1,2384 230 0,2099 1,6833 
1.0 0,9195 1,2933 | 1,0 0,9261 12841 | | 
2,0 0,6970 1.220 | 20 0,6984 1,7184 | 


Zahlentafel 6. Instabiler 


Bereich nullter und erster Ordnung. 


2—5.10-? 


1,58506 | 1,1098 1,1117 | 25 | — 1,4525 | 0,622 1,566 
1,58507 | 1,1088 1,1127 | 21 0,8926 1,275 
1,58510 1,1073 1,1142 |  2,00655 1,0770 1,0799 
158520 | 1,1045 1,1171 2,006578 1,07768 1,07956 
1,58530 | 1.1026 1,110 | 230 — 1,0773 
1,5855 1,0998 11218 15 0,7242 
1,5860 1.0947 1,120 | 10 — 0,4137 
1,5880 1,0835 1.136 | 08 — 0,3096 
1,590 1,0745 1,1480 | 05 — 0,1861 
1,60 1080| 11765 | 08 — 0,1317 
1,70 00390 | 1308 | 0 

1,80 | 13880 | 
2.00 | 0,7805 1,5076 
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Bild 6. Dämpfungskurven in den Iustabilitätsbereichen nullter, halber und erster Ordnung für 
a? = d? — 10-%; 10-3; 10-2; 5. 10-2 und 10-1, 


nauigkeit der Rechnung wurde, da jede weitere Annäherung sehr zeitraubend gewesen wäre, 
dadurch abgeschätzt, daß in die benützten Gleichungen der Wert a=d=0 eingesetzt und 
das Ergebnis mit den genauen Werten von Ince verglichen wurde. Man erhält z. B. für 
das instabile Gebiet erster Ordnung folgende Fehler in °% der Inceschen Werte: 


Abweichungen | Abweichungen 
A, 1, 
1 0,087 
1,5 0,409 °Jo 0,285 
2,0 1,530 ° 0,873 ° 
2,5 4,287 | 2,041 °Jo 


Die Werte der Zahlentafeln 4 bis 6 wurden über die Fehlergenauigkeit (39) berechnet, um 
den Verlauf der Kurven an sich, besonders der flachen Minima, deutlich zu machen. Die 
Fehler sind ‘einseitig, d. h. alle A-Werte sind um die angegebenen Fehler zu weit nach .rechts 
verschoben. 


C. Zusammenfassung und Vergleich der erhaltenen Resultate mit den bisher 
vorliegenden. 


Wir erhalten als wichtigstes Ergebnis, daß die partikulare Lösung für alle jene Para- 
meterwerte A und y der Differentialgleichung eine beschränkte Funktion ist, für die die Resonanz- 


bedingung »+i + =0 (bei Abwesenheit von Dämpfung) bzw. +tu+ ö) =0 


(bei Anwesenheit von Dämpfung) nicht erfüllt ist. Für jene Kurven der Parameterwerte A 
und y, für die die Resonanzbedingung erfüllt ist (Kurven der Resonanzstellen) geht bei Ab- 
wesenheit von Dämpfung die Resonanzlösung im Fall a), wenn also die Störfrequenz weder 
gleich einem ganz- noch halbzahligen Vielfachen der Frequenz der Belegungsfunktion ist, 


mit x gegen Unendlich, während sie in den Fällen £) und y), also für L=5, wenn also die 
Frequenz der Störfunktion gleich der ganzen oder halben Frequenz der Belegungsfunktion 
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bzw. gleich einem ganzzahligen Vielfachen der ganzen oder einem ungeraden ganzzahligen 
Vielfachen der halben Frequenz der Belegungsfunktion ist, mit x? gegen Unendlich geht. Bei 
Anwesenheit von Dämpfung gibt es im Falle a) keine „scharfen“ Resonanzstellen mehr, in 
den Fällen 8) und y) geht die Lösung dann nur mehr mitt x gegen Unendlich. 


Wir haben unsere Untersuchungen der Übersichtlichkeit halber auf die, inhomogene 
Mathieusche Differentialgleichung beschränkt. Die wesentlichen Züge der erhaltenen Re- 
sultate bleiben aber auch für die inhomogene Hillsche Differentialgleichung bestehen. 


Wir wollen noch kurz angeben, was an Ergebnissen. über die inhomogene Hillsche 
(und somit auch Mathieusche) Differentialgleichung bereits in der bis jetzt erschieneneu 
Literatur vorhanden ist. Sehr groß ist,die Merge jener Arbeiten, die sich mit der Hillschen 
Differentialgleichung mit Störfunktion befassen, deren Frequenz gleich der ganzen Frequenz 
bzw. gleich der halben Frequenz der Belegungsfunktion ist. Es handelt sich hier fast immer 
um die Beantwortung technisch wichtiger Fragestellungen. Die benutzten Ansätze sind 
eigentlich immer richtig; es fehlt aber im allgemeinen die genaue Herleitung dieser Ansätze. 
(Siehe auch Erdelyi [6] sowie die dort am Ende der Arbeit zusammengestellte Literatur.) 
Dagegen enthält die Arbeit [6] von Erdelyi eine Herleitung der Form des partikularen 
Integrals einer inhomogenen Hillschen Differentialgleichung mit Dämpfungsglied und einer 
periodischen Störfunktion, deren Frequenz weder ein ganzzahliges der ganzen noch ein un- 
geradzahliges Vielfaches der halben Frequenz der Belegungsfunktion ist; und zwar erfolgt die 
Ableitung der Form des partikularen Integrals gleichfalls über die hier benutzte Methode der 
Variation der Konstanten. Erdelyi erhält genau wie wir das Resultat, daß diese Lösung überall 
in den Gebieten eine beschränkte periodische Funktion zweiter Art ist. Dann aber schließt E. 
weiter aus der Form der Fundamentalintegrale der homogenen Gleichung auf den ursprüng- 
lichen Grenzkurven (Grenzkurven in 2; a=0), daß die Lösung %part hierfür einen mit x 
explizit gebenden Anteil aufweist, also für die Wertepaare A und y der ursprünglichen Grenz- 
kurven unbeschränkt wird, ein Ergebnis, das sich, wie die vorliegende Untersuchung zeigt, 
nicht aufrechterhalten läßt. Die Lösung %»art der inhomogenen Mathieuschen Differential- 
gleichung für diese Störfrequenzen mit (und auch ohne) Dämpfungsglied ist auf diesen Kurven 
beschränkt. Die Resonanzstellen der Lösungen hängen, wie wir nochmal betonen möchten, 


nur von der Erfüllung der Resonanzbedingung (+ u+i+6 + ir)=0 und nicht von der 
Lage der Wertepaare A und y in den Gebieten oder auf den Grenzkurven ab. 


Außer der von E. im Anhang der Arbeit [6] zusammengestellten Literatur über die in- 
homogene Hillsche Differentialgleichung möchten wir noch vor allem auf die Behandlung 
der inhomogenen Hillschen Differentialgleichung, wie sie sich in den Arbeiten über das 
Problem der periodisch veränderlichen Massenträgheit bei Kurbelwellen findet [7], hinweisen. 
(Dieselbe Methode der Behandlung einer inhomogenen Mathieuschen Differentialgleichung 
findet sich auch bei Mettler [8].) Dort wird auf einem anderen Wege als auf dem hier 
eingeschlagenen die Form der Lösung %Ypart für die Fälle 8) und y) abgeleitet und der Satz 
aufgestellt, daß die Resonanzstellen der partikularen Lösung bei jenen Werten der Parameter 
der Differentialgleichung liegen, für die die Lösung der homogenen Differentialgleichung eben- 
falls eine periodische Funktion mit der ganzen oder halben Frequenz der Belegungsfunktion ist. 


Schrifttumverzeichnis. 


[1] K. Klotter: Forsch. Ing.-Wes., Bd. 12 (1941), S. 209. 


[2]M. J. O. Strutt: J,amesche, Mathieusche und verwandte Funktionen in. Physik und Technik. Er- 
gebnisse der Mathematik und ibrer Grenzgebiete. Bd. I, Heft 3, Berlin 1932. Siehe auch das 
Literaturverzeichnis am Ende des Buches. 


[3] E. L. Ince: Proc. Roy. Soc. Edinburgh, Bd. 52 (1931/32), S. 355. 

[4] K. Klotter und G. Kotowski: Z. angew. Math. Mech., Bd. 19 (1939), S. 289. 

[5] E. Mettler: Mitt. Forsch.-Anst. Gutehoffnungshütte-Konzern, Bd. 8 (1940), S. 1. 

[6] A. Erd&lyi: Arch. Elektrotechn., Bd. 29 (1935), S. 473. 

[7] F. Kluge: Ing.-Arch.. Bd. 2 (1931), S. 119. Siehe auch die dort angegebene Literatur. 

[8] E. Mettler: Forschungshefte aus dem Gebiet des Stahlbaues, Heft 4 (1941), 8. 1. 466 


). 
n 
e 
Ss 
| 
zZ, 
0 
)- 
24 
t, 
e 
n 


230 Graf, Affine Transformationen durch Doppel-Photographie zum rn 


Affine 'Transformationen durch Doppel-Photographie. 
Von Ulrich Graf in Danzig. 


Bei der Diapositiv- Herstellung von Schaubildern, Diagrammen, Nomo- 
-grammen, bei der Reproduktion von Zeichnungen aus Büchern u. ä. müßte 
eine in Größe und Format irgendwie vorliegende Urzeichnung auf das 
Diapositiv-Format (beliebiger Größe) umgezeichnet werden, um daraus 
phototechnisch durch ähnliche Verkleinerung das (voll ausgefüllte) Dia- 
positiv zu gewinnen. Unter Vermeidung der Umzeichnung läßt sich rein 
reproduktionstechnisch eine Verkleinerung in verschiedenen Maßstäben 
in Höhe und Breite (also eine affine Transformation) vornehmen, die die 
unmittelbare photomechanische Reproduktion auf beliebige Formate ermög- 
licht. Die bei diesem Verfahren mafßgebenden geometrischen Gesetzmäßig- 
keiten werden untersucht. 


1. Das Photo einer ebenen Figur liefert, sofern von dem trivialen Fall der ähnlichen 
Reproduktion abgesehen wird, ein perspektives Bild der Urzeichnung, bei dem die ur- 
sprünglich parallelen Geraden zu einem Fluchtpunkt zusammenlaufen, die bekannten per- 
spektiven Verkürzungen auftreten usw. Durch ein erneutes Photographieren dieses Photos 
kann unter besonderen Bedingungen erreicht werden, daß in dem zweiten Photo die perspek 
tiven Eigenschaften‘ kompensiert werden und durch die iterierte Photographie ein affin 
transformiertes Bild der Ausgangszeiehnung gewonnen wird, bei dem also die Parallelität und 
die Teilverhältnisse erhalten bleiben. Dieses „fluchtpunktlose Doppel-Photographieren“ wird 
im folgenden geschildert; es ist eine einfache Anwendung des Satzes: Sind zwei Bilder per- 
spektiv zu einem dritten mit derselben Geraden als Verschwindungsgeraden, so sind sie unter- 
einander affin. Die benutzten Photographien sind stets als exakte Zentralperspektiven an- 
genommen. 

2. Die erste Photographie der Urzeichnung {P,} erfolge unter einem Neigungswinkel a, 
der Plattenebene gegen die Zeichenebene (Bild 1, Achsenschnitt). Alle Geraden in {P,}, die 
sich auf der Verschwindungsgeraden v, (Hauptverschwindungspunkt V,) schneiden, werden 
im Bilde {P} parallel, und alle Geraden in {P,}, die parallel sind, haben Bildgeraden in {PB}, 
die sich auf der Fluchtgeraden g, (Hauptfluchtpunkt F) schneiden. Ist d, die Objektiv- 
Entfernung von der Zeichenebene und f die Bildweite der Kammer, sind ferner H, und 


die Hauptpunkte von {P,} und {‘B}, so gilt für den Abstand H, V, der Verschwindungsgeraden 
vom Hauptpunkt H, in {P,} und den Abstand 5 F der Fluchtgeraden vom Hauptpunkt 9 


in {P} 


MV,=d,:ctga, und ÖF=f:ctga,. 


Beim Einstellen der Scharfabbildung ist die Scheimpflug-Bedingung zu beachten, der zufolge 
die Ebene des Objektivs 0, durch die Spurgerade s, (Hauptspurpunkt S,) der Plattenebene 
mit der Zeichenebene laufen muß. 

3. Ist das Zwischenphoto {B} der Urzeichnung {P,} in der geschilderten Weise plıoto- 
technisch hergestellt, so läßt sich aus ihm durch Photographieren in genauer Umkehrung des 
unter 2. geschilderten Ganges die Urzeichnung {P,} wiedergewinnen. Diese Aufnahmelage 
zeigt Bild 2, die sich von Bild 1 nur durch die Stellung und die Ausziehart unterscheidet. 
Die Bilder 1 und 2 können — abgesehen von dem skizzenhaft angedeuteten Photogerät — 
miteinander zur Deckung gebracht werden. 

Wird dagegen das Zwischenphoto {PB} von irgendeiner anderen Stelle O aus unter einem 
Neigungswinkel a aufgenommen und dabei die Forderung erfüllt, daß die Verschwindungs- 
ebene dieser neuen Aufnalıme die Fluchtgerade g, von {P} enthält (d. h., daß die Platten- 
ebene der neuen Aufnahme parallel ist zu der Ebene durch g, und das Objektivzentrum 0), 
so entsteht als Photo des Photos {P} ein neues Bild {P}, Bild 3. Wiederum sind dabei die 
Bilder aller Geraden in {P}, die sich auf g, schneiden, in {P} parallel, oder, anders aus- 
gedrückt, in den beiden Bildern {P,} und {P} geht die unendlich ferne Gerade perspektiv 
aus ein und derselben Geraden g,„ eines dritten Bildes {®} hervor. Dann sind die beiden 
Bilder {P,} und {P} zueinander affın, und da {P,} die vorgelegte Urzeichnung war, ist von 
ihr in {P} photomechanisch eine Affinfigur gewonnen. Bild 4 zeigt zusammen die beiden 
Photolagen {P,} und {P}, durch die aus dem nach Bild 1 erzeugten Zwischenphoto {P} so- 
wohl die Originalzeichnung {P,} als auch ihr affines Bild {P} erzeugbar ist. Wird an Stelle 
des Öriginalphotos {P} eine A-fache Vergrößerung benutzt, so ist überall f durch Af zu 
ersetzen. 

4. Zur analytischen Erfassung dieses geometrischen Sachverhaltes, den Bild 5 in anschau- 
licher Skizze zeigt, wird ein X YZ-Achsenkreuz benutzt, dessen X Y-Fbene das Zwischen- 
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Bild 1. Die Herstellung des Zwischenphotos {®} 


Bild 4. Die Aufnahmelagen für die Urzeichnung 
und das Affinbild. 


Bild 2 (liuks). Zwischen- 
"photo und Urzeichnung. 


N 


Bild 3 (links). Zwischen- 
photo und Affinbild. 


Bild 5 (rechts). Der Zu- 
sammenhang der Auf- 
nahmedaten. 


photo {P} enthält und dessen Y-Achse parallel zu der Verschwindungsgeraden g, im Ab- 
stand f-ctga, verläuft. Die Zentren O, und O der von diesem Zwischenbild herstellbaren 
Photos {P,} und {P} haben die Koordinaten 0, (X=0, Y=0,Z=f) und 0(X=X,,Y=Y,, 
Z=Z,). Die Neigungswinkel der beiden Plattenebenen //, und /] seien «a, und a, die Bild- 
weiten d, und d. Gemäß der Affinitätsforderung müssen die Verschwindungsebenen der beiden 
Aufnahmen die Verschwindungsgerade g,„ enthalten, d. h. es muß 
oder .. . (l) 
sein. In den Bildebenen /J, bzw. II werden Koordinatenkreuze £&,, n, bzw. &, 7 benutzt, 
deren Ursprung der Punkt A, bzw. A ist, der als Bild des Hauptpunktes ö des Zwischen- 
photos {P} entsteht. Die n,- bzw. n-Achse verläuft parallel zur Verschwindungsgeraden 4,, 
die bzw. &-Achse windschief-senkrecht zu ihr. 
Die Bildebene /J hat die Gleichung 

Der Strahl RO von einem Punkt P(X=x, Y=y,Z=0) der X Y-Ebene durch das Zentrum 
O(X,, Y,, Z,) trifft die Plattenebene /J im Bildpunkt P mit den Koordinaten 

d-(X, — x) 
+4 
d-(Y,—y) 
— x)-sina+Z,:-cosa’ 
d.Z, 

— x)-sina+ Z,-cosa 
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Wird als Sonderfall an Stelle des allgemeinen Punktes ® der Hauptpunkt 9 (x=0, y—0) 


in (B} gewählt, so wird man auf den Ursprung A (X,, Y4, Z,) des £n-Systems in der Ebene /] 
geführt. Es ist demnach 


Die Plattenkoordinaten &,n des Punktes P in // sind 


sie ergeben sich nach (1), (2) und (3) zu 
y:f—x-Y,-tga, 
Die Koordinaten £&,, n, des Bildpunktes P, in //,, der als Bild von ® bei der Aufnalıme 


mit dem Zentrum 0, entsteht, erhält man aus (5) für ,=Y,=0, Z,=f, a=«a, d=d, 
E=&,,n=n. Es ergibt sich 


und die Elimination von x und y aus (5) und (6) liefert die gesuchte affine Bindung zwischen 
den Ebenen //, und IT. Die Ausrechnung liefert 


& 


Es verhalten sich 


"sina’sina, 

d. h. wie die Abstände der Zentren O und O, von der gemeinsamen Verschwindungsgeraden g, 
(Bild 4). Der Koeffizient a, ist proportional der Höhenverschiebung Y,. Die in (7) gewonnene 
Sonderform der affinen Bindung zwischen /T, und /I ist durch die geeignete Wahl der Platten- 


koordinaten erreicht; eine zusätzliche Verschiebung und Drehung einer Platte in ihrer Ebene 
führt auf die allgemeinste Form der Affinität. 


5. Die allgemeinen Formeln (7) der affinen Doppel-Photographie lassen folgende praktisch 
wesentlichen Sonderfälle erkennen. 


a) Im Sonderfall «=0 und damit auch a=0 ergeben sich ähnliche bzw. kongruente 
Abbildungen. In der Tat ist 
im _ A 
a. >oSina Z, 
also 


b) Ist ,=+0, so wird die Affinität zur Ähnlichkeit für a, —=b,, a,=0. Das tritt’ ein für 


sina sina, 


Der geometrische Ort der Objektiv-Zentren O ist in diesem Falle der Kreis um F 
durch O, (Bild 4')). Bei jeder solchen Aufnahmelage ist wie eben 


1) Vgl. U. Graf: Uber das Photo eines Photos, Z. angew. Math. Mech. Bd. 21 (1941), S. 183 bis 189. 
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(3) 


(4); 


(7). 


für 
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e) Für Y,=0, a=a, ist das Aufnahmegerät ohne Höhenänderung in seiner Verschwin- 
dumgsebene verschiebbar; die Balın von O ist die Gerade O0, F (Bild 4). In diesem 
Falle ist 


Wird das so entstandene Bild noch im Maßstab @: d, vergrößert, so hat man in der 
Transformation 


denjenigen Sonderfall der Affinität vor sich, der, durch den Parallelschatten realisiert. 
in der Darstellenden Geometrie eine besondere Rolle spielt ?). 


6. Bei der praktischen Durchführung des geschilderten affinen Doppel-Photo-Verfahrens 
könnte der Umstand erschwerend wirken, daß das auf der Plattenebene // erzeugte Bild weit 
von dem Hauptpunkt H in II entfernt liegt. ‚Erstrebenswert ist es, das Bild A des Haupt- 
punktes ö im Zwischenphoto {PB} in die Nähe des Hauptpunktes H der Plattenebene /T fallen 
zu lassen. 

Die Raumkoordinaten Xy, Ya, Zu von H sind 


Xy=X,+d-sina, Zu= 2, + d:cosa, 


und nach (1) und (4) ergeben sich die Plattenkoordinaten von H zu 


X,:cosa—Z,-sina , X. (fr X,)—Z 

(8). 
d 


Die Forderung &4 = const liefert als geometrischen Ort für die Objektivzentren O(X,, Y,, Z,) 
den Zylinder (Bild 6) 


oder 


der die Y,-Achse und die Verschwindungsgerade g, als Mantellinien enthält und dessen Grund- 


kreis in der X, Z, Ebene über der Strecke 9 F einen Bogen mit dem Peripheriewinkel a+®P 
aufspannt, sofern 


p= are 


der auf die Bildweite d bezogene zu &y gehörige Abweichwinkel ist. 


Die Forderung 7,5 = const liefert als geometrischen Ort für die Objektivzentren O(X,, Y,,Z,) 
ein gerades Konoid mit der Verschwindungsgeraden g, als Achse (Bild 7). Die gerad- 
linigen Erzeugenden stehen senkrecht auf g, und durchstoßen den Zylinder um -die Achse g, 
mit $F=f.cetga, als Halbmesser in einer Ellipse, die aus dem Zylinder von der Ebene 


Y=tgy-Z, 
ausgeschnitten wird, sofern 


y g d 


wieder der auf die Bildweite d bezogene zu 7 gehörige Abweichwinkel ist. Die übliche 
Parameterdarstellung dieses Konoids lautet 


Z=u-sinv . . . (10), 


wobei » der Winkel der -Erzeugenden gegen die X, Y,Ebene und « die Länge der Er- 
zeugenden sind. 


2) Ein Verfahren, Parallelprojektionen mit einem Photogerät durch gewisse Abblendungen zu erzeugen, be- 
schreibt E. Gotthardt: Zur Frage der Projektionszentren in der Bildmessung ; Bildmessg. u. Luftbildwes. Bd. 16 (1941), 
S. 97/98. Damit können auch räumliche Gegehstände parallelprojiziert abgebildet werden, jedoch sind die Diuggrößen 
auf einen zylinderförmigen Raum vom Durchmesser des Objektives beschränkt. 
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Bild 6. Der Zylinder 
—=d.tgp=const. 


Bild 7 (rechts). Das Konoid 
tg v= const. 


20 


Die vollständige Fläche, die sich für |77|=const durch rückwärtige Verlängerung der 
Erzeugenden (also durch zusätzliche Spiegelung des Konoids in Bild 7 an der X, Y,-Ebene) 
ergibt, ist eine Regelfläche vierter Ordnung mit der Gleichung 


(X, —f- ... 0.0.0. MM. 


Nach der Einteilung von R. Sturm) ist sie von 7. Art; sie besitzt die doppelt-zählende 
Verschwindungsgerade g, und die doppelt-zählende unendlich ferne Gerade der X, Z,- Ebene 


als Leitgeraden und in der X,-Achse eine Doppel-Erzeu- 
gende. Jede Ebene durch die X,-Achse schneidet die 
Fläche außerdem in einer Ellipse, getragen von einem 
Zylinder um g, als Achse, dessen Halbmesser proportional 
zu dem Anstieg der Schnittebene gegen die X, Y,- Ebene 
ist; der Proportionalitätsfaktor ist f- etga,-tgw. 

Als geometrischer Ort der Öbjektivzentren O(X,, Y,,Z,) 
für alle Aufnahmen, bei denen die Hauptpunktkoordinaten 
£4, nu fest bleiben, ergibt sich die Schnittkurve des 
Zylinders (9) mit dem Konoid (10), die in Bild 8 für 
Y,=0 gemäß dieser Erzeugung konstruiert ist. (Für die 
affine Doppel: Photographie ist nur der Kurventeil Z, > 0 
von Bedeutung.) Die Parameterdarstellung der Kurve 
lautet’nach (9) und (10) 


cos v 
) 
Y,=f-etga,-tgy-sinv, 


sin v 


wobei » die oben erklärte Bedeutung hat. Diese Kurve 
vierter Ordnung ist der wesentliche Bestandteil der Ge- 
samtdurchdringungskurve achter Ordnung der Regelfläche 


20 


Bild 8. Die Raumkurve Sg > eonst, 
N > eonst als Durchdringung von Zylinder 
und Konoid. 


(11) und des Zylinders (9), die außerdem noch die doppeltzählende Gerade g, und das kon- 
Jugiert-imaginäre Geradenpaar X5+Z;=0 der unendlich fernen Ebene enthält. 


3) R. Sturm: Liniengeometrie I, Leipzig 1892, 8. 52 bis 61. 
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Die Kurve (12) läßt sich als der Schnitt des Zylinders (9) mit dem Rotations 
ellipsoid 
) 


Ku=f-cga, 


erzeugen und erweist sich so als eine aus der konstruktiven Technik geläufige Durchdringungs- 
kurve. Die Rotationsachse ist die Zylindermantellinie, die der Y,Achse diametral gegenüber- 
liegt. Ellipsoid und Zylinder berühren einander im Punkte $, einem Doppelpunkt der Durch- 


dringungskurve. Die Scheitel des Ellipsoids liegen in Höhe der höchsten und tiefsten Er- 
zeugenden des Konoids (10). 


Das Büschel von Flächen zweiter Ordnung durch die Kurve (12) enthält außer dem 
Kreiszylinder (9) noch den parabolischen Zylinder 


dessen Scheitelmantellinie die Tangente in ö an den Grundkreis des Zylinders (9) und dessen 
Parameter 


ist. Die Raumkurve (12) läßt sich daher als Durchdringung der Zylinder (9) und (14) gewinnen 


(Bild 9). 


Bild 9. Die Raumkurve $,,—=const, 7,, =const als Durch- Bild 10. Der Raumteil |$,,|=const, |n,,|=eonst. 
dringung von geradem Kreiszylinder und parabolischem 
Zylinder. 


Als Ergebnis der Überlegung folgt: 


Soll bei dem affinen Doppelphoto-Verfahren das Bild A des Hauptpunktes 5 im Zwischen- 
photo (P} auf der Plattenebene innerhalb eines Rechtecks mit den Seiten und 
um den Plattenhauptpunkt H als Mittelpunkt liegen, so muß das Objektivzentrum innerhalb 
des in Bild 10 wiedergegebenen Raumteiles liegen. Er wird begrenzt von einem Zylinder- 
BEP über SF als Sehne, wobei die bestimmenden Peripheriewinkel der Kreisbogen 


3 +arctg Eu sind, und zwei Konoiden (10) mit Grenzerzeugenden in den Höhen ,=+9HF° 1 


d' 
Die Kurvenkanten sind Teile der geschilderten Raumkurven vierter Ordnung. 
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Die phototechnisch zu- 
teiles erweist die praktische 
- - - .. . 
Durchführbarkeit des geschil- 
+ 1 

derten Doppelphoto - Ver- 

N 2 fahrens. 

| | das Bild von in den Haupt- 
punkt H selbst. In diesen 

3-2 Falle entartet der Raumteil 

Ho zu dem Kreis 

der F als Durchmesser be- 
7 sitzt, ein Sachverhalt, der ele- 
4 mentargeometrisch unmittel- 

bar einleuchtet (Bild 4). Bei 
der Einhaltung dieser beson- 
die allgemeine 7 
Bild 11. Die Urzeichnung (Seitenverhältnis 9: 14). u 8 n Formeln ( 
die Gestalt an: 
| d sine, d sin “) 
d, sina d, sina)' 
m 
n 
4 7. Bild 11 zeigt eine Urzeich- 
—u — nach (15) in ein Quadrat (Bild 12) 
4 umgeformt wurde. t 
1-4 war gemäß der Forderung 
der zweite Neigungswinkel a = 33°. 
1%  F%. fahren besonders bequem für 
34-0 ON Et? a=a,, da man in diesem Falle 
. r . 
I mit nur einer Kamıner-Einstellung 
Hu und Kurvenbilder, die nicht wie me- 
4 - . 
an feste Maßstäbe gebunden sind, 

hinaus kann man durch iteriertes 

15 LH} Photographieren nach (7) jede 
4 HH 4H ebene affine Transformation er- 

M zeugen. 475 
2 2 ) Vgl. Handbuch der modernen Re 


Bild 12. Die Affinreproduktion (Seitenverhältnis 1:1). 


1936, 


produktionsteehnik, Frankfurt a. M. 
Bd. I, S. 141. 
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KLEINE MITTEILUNGEN 


Graphische Lösung von Randwertpro- 
blemen bei gewöhnlichen linearen Diffe- 
rentialgleichungen 2. Ordnung. 


1. Fragestellung. 

Das im folgenden beschriebene Verfahren ist ein 
Näherungsverfahren; es wird der Differential- 
gleichung (1) eine Differenzengleichung (4) gegen- 
übergestellt und das zugehörige Differenzen-Rand- 
wertproblem graphisch gelöst !). 

Das vorgelegte Randwertproblem habe die Diffe- 
rentialgleichung 


und die Randbedingungen 


An Stelle der zweiten Randbedingung kann auch 
y (b)=y treten. Über kompliziertere Randbedin- 
gungen vgl. Nr. 4. 

Es seien a, ß, y gegebene Konstanten, a, b end- 
liche Abszissen mit a<b, und r (x). f(x) 
gegebene stetige Funktionen von 7. Striche bei Y 
bedeuten Ableitungen nach r. 

Das Randwertproblem wird angenähert mit dem 
Differenzenverfahren gelöst: bei Teilung des Inter- 
valls 5 in n gleiche Teile der Länge = = 


Y(a)=a, 


mit den Teilungspunkten 


werden Funktionswerte an der Stelle x, durch 
Anhängen des Index bezeichnet, z. B. 9: = (ax); 
ferner bedeutet %7, einen Näherungswert für den 
Wert der Lösungsfunktion Y(z) an der Stelle z;.. 
Der Differentialgleichung wird die Differenzen- 
gleichung 


h? +ri | 
. r;h 
oder mit der Abkürzung 9, =, 


gegenübergestellt. Als Randbedingungen kommen 
MmMuh..... 0) 


hinzu. Die Gl. (4) (5) bilden das zu lösende Diffe- 
renzen-Randwertproblem, nämlich (rn —+ 1) lineare 
Gleichungen für ebenso viele Unbekannte %,, Y1, 


2. Die Sehnenfestpunkte. 


Man könnte das Randwertproblem (4), (5) lösen, 
indem man aus der ersten Gl. (4) %, durch %, aus- 
drückt (4, ist ja bekannt), dann aus der zweiten 
Gl. (4) durch y, ausdrück&usw. Hat man auch 
Y„n als lineare Funktion von %, dargestellt, so 
liefert die Randbedingung für %, den Wert von %,.. 


1) Der Aufsatz ist unter Betonung geometrischer Ge- 
sichtspunkte eine leichte Verallgemeinerung einer Arbeit 
von R. Kappus: „Graphische Lösung der Differential- 
bei beliebigen Anfangs- und 

andbedingungen“, Z. VDI Pd. 87 (1943), S. 26 bis 28, in 
welcher über ein Verfahren von J. Drymael: J. roy. 
aeron. Soc, Bd. 45 (1941) Nr. 362, S. 51 bis 66, referiert wird. 


Es erscheinen so alle 4; als lineare Funktionen 
von %ı 


Y%k Yıt dk (k=1,2,...,n) 


mit a;, 5b; als Konstanten. Denkt man sich bei 


-festem , =a für mehrere verschiedene Werte von 


y, nach (4) die jeweils zugehörigen Werte von %, 
Y3,...,%% berechnet und dargestellt, Bild 1, so 


werden auf den zu 2,,- ..2, gehörigen Ordinaten 
ähnliche Punktreihen ausßeschnitten, d.h. aber die 
Verbindungsgeraden entsprechender Teilpunkte 
schneiden sich in einem Punkt. Insbesondere 
schneiden sich die Verbindungsgeraden der Punkte 
(£j-,%i-ı) und (z,,4;) in einem Punkt, der als 
Sehnenfestpunkt S; bezeichnet werde. ?) 

Ferner werden noch die Tangentenfestpunkte 7; 
eingeführt. Als „Tangente“ im Punkt z,,%; wird 
die Parallele zur “Werbindungslinie der Punkte 
(&;-»%:-ı) und (7j+1,%i+,) definiert. 7; ist dann 
der Schnittpunkt aller Tangenten in z;,y;. Greift 
man unter den zu den verschiedenen %;- Werten 
gehörigen Linienzügen denjenigen heraus, für den 
der Punkt (z;,%;) auf der Verbindungsgeraden von 
S; und $;,, liegt, so stellt diese Verbindungs- 
gerade zugleich die Tangente im Punkt (z;.%;) 
dar, d.h. S;.S;,,.7; liegen auf einer Geraden. 
Bild 2. 


y 5 
l 
Bild 2. 


3.Konstruktion der Sehnenfestpunkte. 


Der erste Sehnenfestpunkt S, ist der gegebene 
Punkt (2,,%0). Wir geben an, wie man bei be- 
kanntem Festpunkt S; den nächsten Festpunkt 
S;+, konstruieren kann. Damit ist nämlich das 
Randwertproblem (4), (5) leicht lösbar: Hat man 
alle Sehnenfestpunkte S,, S,,...,S„ konstruiert, so 
kennt man die letzte Sehne als Verbindungslinie 
der Punkte S„ und (2,,%„) und kann dann die 
vorletzte Sehne und so von rückwärts den ganzen 
Sehnenzug zeichnen. Vgl. in Bild 1, den strich- 
punktierten Linienzug. 

Es seien &;,n; die Koordinaten des Sehnenfest- 
punktes $S;, entsprechend &; + ,,Ni+, die von 


?) Festpunkte werden auch verwendet zur graphischen 
Behandlung der Clape yronschen Dreimomentengleichun- 
gen, vgl. P. Funk, Die linearen Differenzengleichungen 
und ihre Anwendung in der Theorie der Baukonstruk- 
tionen, Berlin 1920, S. 44. 
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Bild 2; die drei Punkte (z;_,,4i-,) (2i,4i) und S; 
liegen auf einer Geraden, d.h. Bild 2. 
— Yi-ı (6) 
Genau so folgt: 
YıMizsı 
...00. 
Yrıdi Mm. 
Mit Hilfe von (6) und (7) kann man y;_, und 


Yi+. aus der Differenzengleichung (4) eliminieren 
und man erhält 


—fir 
Diese Gleichung soll nun für beliebige Werte 
von %; gelten, d.h. es muß der Koeffizient von %; 
und das absolute Glied, jedes für sich ver- 
schwinden: 


1+60i 10 


Bei bekanntem Sehnenfestpunkt S; kann man 
hiernach den nächsten Festpunkt S;,, berechnen, 


nämlich aus (9) die Koordinate &;,, und dann 


aus (10) n;+.. Die Gleichungen lassen eine ein- 
fache graphische Konstruktion, Bild 3, zu. Außer 


gh? 
1 
| | 
Sie £ Ti 
Bild 3. 


der z-y-Ebene‘ verwendet man eine Hilfsfigur, in 
welcher über der z-Achse die Kurven 1—o,1-+0 
und 1+o— „ph? aufgetragen sind, welche bei 
z=:t; die Punkte B,C,D auschneiden. Vom 
Punkte A(z=$;) ausgehend, zieht man zu AB die 
Parallele durch D, die die Gerade z=z;,, inE 
schneidet. Die Gerade CE liefert dann in dem 
Schnittpunkt F mit der z-Achse die gesuchte 
Abszisse &;+.,. Die Richtigkeit der Konstruktion 
rechnet man unmittelbar nach; der ÖOrdinaten- 
unterschied zwischen den Punkten D und E ist 


a, der zwischen den Punkten C und E 


also - und dieser soll nach (9) 
gleich sein. Für die geometrische 


Konstruktion der noch fehlenden Ordinate n;,, 
empfieblt es sich, Bild 4, durch Ziehen der Par- 
allele CG zur Geraden AB eine Abszisse £; mit 
_ 
1-0 1+60 


zu konstruieren. Dann lassen sich (9) und (10) 
in der Gestalt schreiben 


1 h 


Division dieser beiden Gleichungen liefert 


es Eirı 


Bild 4. 


Zeichnet man sich in der z-y-Ebene die Kurve für 
L, die bei z=z; durch den Punkt H geht, Bild 4, 


und bezeichnet man den Punkt mit den Koor- 
dinaten £{;, n; mit J, so liegen nach Gl. (11) die 
Punkte H, J und der gesuchte Sehnenfestpunkt 
S;+, auf einer Geraden. Dieser ist also jetzt 
völlig festgelegt. Der Gang der Konstruktion ist 
durch Pfeile angedeutet. Im Falle r («)=0, also 
0;==0 vereinfacht sich die Konstruktion und man 
erhält eine Figur wie bei R. Kappus a. a. O0. 


4. Andere Randbedingungen. 


Ist an Stelle von y(d)=3 die Randbedingung 

(d)=y vorgegeben, so braucht man den letzten 
Tangentenfestpunkt 7,„, um durch ihn die letzte 
Strecke des Polygonzuges für y mit. der vorge- 
ze Steigung zeichnen und dann die übrigen 

trecken des Lösungszuges wie früher von rück- 
wärts her ermitteln zu können. Wir beschreiben 
allgemein die Konstruktion von 7;, wenn S; und 
+, bekannt sind. 7; habe die Koordinaten £}, n}. 
Da nach Nr. 2 8;,S;,,,7; auf einer Geraden liegen, 
genügt es, die Abszisse &} zu konstruieren. 

Es soll nach der Definition des Tangentenfest- 
punktes gelten 


_Yirı—di-ı 


Jetzt wird eine Rechnung gleicher Art durchge- 
führt, wie sie zu (8) führte. Aus (6) und (12) 
werden die Differenzen —Y;_, und 
berechnet und in die Differenzengleichung (4) ein- 
gesetzt. Es entsteht eine in %; lineare Gleichung, 
die in ; identisch gelten soll, es müssen also 
Absolutglied und Koeffizient von y; einzeln ver- 
schwinden; der Koeffizient von y; gleich Null ge- 
setzt liefert 


Diese Gleichung &®ann nach Bild 5 geometrisch 
edeutet werden. Zeichnet man über der z-Achse 
ie Kurven für 1+0,1,1—g#?, die bei r=z; 

die Punkte A,B,C ausschneiden, so geht man vom 
Punkt D mit der Abszisse =}; aus. Die Par- 
allele zu B D durch den Punkt C trifft die Ge- 
radez=r;,,inE, und die Parallele zu BE durch 
A schneidet im Punkte F die gesuchte Abszisse & 
aus. 
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Bild 5. 


Bei komplizierteren Randbedingungen, z. B. 
y(b)=ö-y(b), kann man zum Ziel kommen, indem 
man die Konstruktion zweimal durchführt, also 
zwei linear unabhängige Lösungen der Differential- 
gleichung konstruiert und aus ihnen durch Linear- 
kombination eine die Randbedingungen erfüllende 
Lösung ermittelt. 


5. Beispiel. 


‚Bild 6 zeigt die Durchführung eines einfachen 
analytischen Beispiels mit bekannter exakter 
Lösung. Das Randwertproblem 


-(143)» 


+ 
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Bild 6. 


hat die Lösung y=1—r°. Bei der groben 
Maschenweite 1-4 hat man vom gegebenen An- 


fangspunkt S, (r=-—1;y=0) ausgehend die drei 
weiteren Sehnenfestpunkte S,, S,, S, zu konstruieren 
und erhält in dem strichpunktierten Linienzug die 
Näherungslösung des Differenzenverfahrens. 


Karlsruhe. L. Collatz. 473 


BUCHBESPRECHUNGEN 


Dr. HORST VON SANDEN, o. Prof. a. d. Techn. 
Hochschule Hannover, Praxis der Diffe- 
rentialgleichungen, eine Einfüh- 
rung. 100 S. m. 20 Abb. Berlin 1943, Verlag 
Walter de Gruyter & Co. Preis brosch. 5 M. 


Das vorliegende Buch, das den Inhalt von Vor- 


lesungen wiedergibt, die Verf. seit längerer Zeit 
an der Techn. Hochschule Hannover hält, behan- 


delt in der Hauptsache Näherungsmethoden zur 
Lösung von Differentialgleichungen. Der erste 
Teil beschäftigt sich mit Anfangswertproblemen 
Differentialgleichungen. Als Metho- 
en zur Lösung von Gleichungen erster Ordnung 
findet man das Isoklinenverfahren, die beiden 
Adamsschen Verfahren, die Iteration und das Ver- 
fahren von Runge-Kutta. Die erforderlichen Ihter- 
polationsformeln, wie die Formeln für die letzte 
Methode werden nicht abgeleitet, sondern betreffs 
der Herleitung wird auf andere Bücher verwiesen. 
Weiter werden für Gleichungen zweiter Ordnung 
und fürSysteme von Gleichungen die zeichnerische 
Integration und das Adamssche Verfahren behan- 
delt, wie die Lösung durch Potenzreihen. 

Im zweiten Teil beschäftigt sich Verf. mit Rand- 
wertproblemen bei Differentialgleichungen der 
Form u” +Ap(a)y=q(re). Zunächst wird an der 
Gleichung mit konstanten Koeffizienten die Frage- 
stellung erläutert, weiter auf Gleichungen mit 
variablen Koeffizienten der Ansatz von Rayleigh 
und die Iterationsmethode angewandt, kurz auf 
die Greensche Funktion eingegangen und schließ- 
lich die Methode von Ritz auseinandergesetzt. 


Nicht behandelt wird also das Differenzenverfahren, 
das gelegentlich dann von Nutzen sein kann, wenn 
man schnell einen Überblick über den Verlauf der 
Lösung gewinnen will. An einer beträchtlichen 
Anzahl von Beispielen wird die Anwendung der 
einzelnen Methoden erläutert. Die Darstellung ist 
sehr klar und ziemlich ausführlich; ihr kommt 
überall die reiche Unterrichtserfahrung des Verf. 
zugute. 


Auf einige Einzelheiten sei kurz hingewiesen, 
z. B. auf die Ausführungen über „exakte“ und 
„Näherungslösungen“, wo man hinzufügen möchte, 
daß diese nur dann als gleichwertig angesehen 
werden sollten, wenn man den Fehler der Nähe- 
rungen genügend genau abschätzen kann, was lei- 
der oft nicht der Fall ist. — Das Problem der 
Knickung eines Stabes unter seinem Eigengewicht 
löst Verf. durch Potenzreihen und sagt dann, daß 
für diese „Symbole noch nieht eingeführt und 
Tabellen noch nicht vorhanden (?)“ sind. Die 
Lösung läßt sich durch Besselsche Funktionen 
darstellen, und zum mindesten die Nullstellen kann 
man aus Tabellen, z. B. im Jahnke-Emde, entneh- 
men. — Zu dem mehrfach mit Näherungsmethoden 
behandelten Beispiel der Schwingung einer Saite 
mit abschnittsweise konstanter Massenbelegung _ 
wäre wohl zu erwähnen, daß sich die Eigenwerte 
aus einer einfachen transzendenten Gleichung er- 
gr, so daß man diese leicht eingabeln kann. 

ur Ableitung dieser Gleichung ist an Mathematik 
nur das erforderlich, was Verf. im letzten Ab- 
schnitt seines Buches entwickelt. — Das auf S. 70 
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angegebene Verfahren zur Berechnung höherer 
Eigenwerte wird an keinem Beispiel erläutert. Ob 

o praktisch brauchbar ist, erscheint mir zweifel- 
aft. 

Diese Bemerkungen beziehen sich auf neben- 
sächliche Einzelheiten. Als Ganzes ist das Buch 
vortrefflich geeignet, Interesse für die rechnerische 
Behandlung von Differentialgleichungen zu wecken 
und in dieses Teilgebiet der praktischen Analysis 
einzuführen. 


Dresden. 


Dr. GERHARD KOWALEWSKI, Prof. a. d. Deut- 
schen Karls-Universität Prag, Einführung in 
die Determinantentheorie einschließ- 
lich der Fredholmschen Determinanten. 
3., verb. u. erweit. Aufl. 320 S. Berlin 1942, Ver- 
lag Walter de Gruyter & Co. Preis geb. 15 M. 


Die Determinante ist von Leibniz und Cramer 
ersonnen worden zur Auflösung von linearen 
Gleichungen. Es hat sich zwar herausgestellt, daß 
.man die Bedeutung der Determinanten für diese 
Frage zunächst überschätzt hat, sowohl was ihre 
begriffliche Klärung als auch was ihre numerische 
Beherrschung betrifft; dafür aber haben sich die 
Determinanten an zahlreichen anderen Stellen als 
nicht mehr entbehrliches Hilfsmittel bewährt. 
Gerade von dieser Vielfalt ihrer Verwendungs- 
möglichkeit gibt das bekannte Buch, das jetzt in 
dritter Auflage vorliegt, einen glänzenden Ein- 
druck. Es wendet sich an Leser, denen eine Ein- 
führung, die nur das Nötigste bringt, nicht genügt. 
Dementsprechend findet man in der zweiten Hälfte 
des Buches Gegenstände wie: Orthogonale Deter- 
minanten, Resultanten und Diskriminanten, lineare 
und quadratische Formen, Funktionaldeterminan- 
ten, Wronskische und Gramsche Determinanten, 
lineare Integralgleichungen, Elementarteilertheorie. 

Das Werk wird gewiß auch in der neuen Auflage 
seine alte Anziehungskraft zeigen. 


Dresden. 


Dr. K. SCHWIDEFSKY, wissenschaftl. Mitarbei- 
ter d. Opt. Werke Zeiß, Jena, Einführungin 
die Luft- und Erdbildmessung. 3. er- 
weiterte u, verbess. Aufl. VI + 176 S. m. 85 Abb., 
3 schwarzen und 2 farbigen Tafeln im Text, 
1 schwarzen Tafel, 1 farbigen Brille und 2 Stereo- 
bildern im Anhang. Leipzig und Berlin 1942, Ver- 
lag B. G. Teubner. Preis geb. 7,80 M. 

Der Fortschritt in der Entwicklung der Photo- 
grammetrie, diesem neuen besonders für topogra- 
phische Zwecke wichtigen Meßverfahren ist noch 
immer recht lebhaft, wie sich aus der bereits nach 
2 Jahren erscheinenden weiteren 3. "Auflage des 
Buches von Schwidefsky „Einführung in die Luft- 
und Erdbildmessung“ ergibt. (Siehe Z. angew. 
Math. Mech. Bd. 18 (1938) S. 259 und Bd. 20 (1940) 
S. 240. Die neue Auflage hatte zwar weniger in- 
strumentelle Neuerungen zu berücksichtigen, als 
solche auf dem Gebiet der mathematischen Grund- 
lagen und der praktischen Anwendung. Die Haupt- 
aufgabe der Photogrammetrie oder die räumliche 
Doppelpunkteinschaltung wird neu für Schrägauf- 
nahmen behandelt, die gefährlichen Örter bei der 
Hauptaufgabe sind erwähnt, der Abschnitt über 
innere Orientierung von Meßkammern und die 
Fehlertheorie der Aerotriangulation sind vertieft, 
ferner wurden Abschnitte über Fehlerfortpflanzung 
in gleichförmigen gestreckten Dreiecksketten und 
die Drehzerrung (günstigste Netzeinpassung) be- 
rücksichtigt. Nach der praktischen Seite neu ist 
vor allem der Abschnitt über die Nahbildmessung 
und das Beispiel zur Aerotriangulation. Durch 


Willers. 468 


Rellich. 443 


diese wichtigen Neuerungen hat das _ treffliche 
Buch an Wert noch gewonnen. 


Hannover. Finsterwalder. 457 


Dr.-Ing. P. WERKMEISTER, o. Prof. em. a. d. 
Techn. Hochschule Dresden, Vermessungs- 
kunde IH, Trigonometrische und ba- 
rometrische Höhenmessung, Tachv- 
metrieund Topographie. 4. Aufl. (Samm- 
lung Göschen Bd. 862). 147 S. m. 64 Fig. Berlin 
1942, Verlag Walter de Gruyter & Co. Preis 1.62 M. 

Es ist recht erfreulich, daß es der Verlag de 
Gruyter &' Co. ermöglichen konnte, trotz aller 
kriegsbedingten Schwierigkeiten nach dem I. Band 
(Stückmessung und Nivellieren, 7. Auflage) und 
der Fortsetzung in Band II (Messung von Horizon- 
talwinkeln, Festlegung von Punkten im Koordi- 
natensystem, Absteckungen, 4. Auflage) nunmehr 
auch den abschließenden III. Teil (Trironome- 
trische und Barometrische Höhenmessung, Tachy- 
metrie und Topographie) der Vermessunzskunde 
von P. Werkmeister als Band 862 der Sammlung 
Göschen in 4. Auflage herauszugeben, und zwar 
in bekannt guter Drucklegung und Ausstattung, so 
daß den bewährten vorzüzlichen Ausführungen des 
Verfassers trotz der Kriegszeit der verdiente 
äußere Rahmen in ansprechender Form gesichert 
ist. 

Freiberg-Sa. H. Müller. 456 

Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten}: 


Dr. A. TIMPE, o. Prof. a. d. Techn, Hochschule 
Berlin, Höhere Mathematik I, Tl. Auto- 
graphierte Vorlesungen. 84 u. 130 S. Stuttgart 
und Berlin 1942, Verlag W. Kohlhammer. Preis 
brosch. Teil 13 M, Teil II 3,50 M. 


Dr. ROBERT SAUER, o. Prof. a. d. Techn. 
Hochschule Aachen, Theoretische Ein- 
führung in die Gasdvynamik. VII+ 
146 S. mi. 99 Abb. Berlin 1943, Springer-Verlag. 
Preis brosch. 12 M. 


GREGOR WENTZEL, Prof. a. d. Univ. Zürich, 
Einführung in die Quantentheorie 
der Wellenfelder. IV+209 S. Wien 1943, 
Verlag Franz Deuticke. Preis brosch. 20 M. 


P. RAETHJEN, Einführungindie Phy- 
sikderAtmosphäre,Bd.Il, Meteorolo- 
gische Aerodynamik. (Hamburger Mathe- 
matische Einzelschriften, 33. Heft/1942.) VII + 
257 S. m. 57 Abb. Leipzig und Berlin 1942, Ver- 
lag B. G. Teubner. Preis geb. 10 M. 
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